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Numerical Calculation of Optimum a for a 
Germania-Doped Silica Lightguide 


By G.E. PETERSON, A. CARNEVALE, U. C. PAEK, andJ. W. FLEMING 
(Manuscript received September 11, 1980) 


Using an exact numerical solution to Maxwell’s equations, we 
determine the optimum refractive index profile parameter a for a 
germania-doped silica lightguide. Our results agree closely with the 
earlier work of Olshansky and Keck and with the work of Marcatilt. 
However, differences exist that may be important in the manufacture 
of very high bandwidth lightguides. 


!. INTRODUCTION 

It is now well established that a lightguide with a near parabolic 
index profile can have single frequency bandwidths of 1 GHz-km or 
greater.’ This is expected to happen if the profile is smooth, the 
central dip is narrow or small, the profile parameter is properly chosen, 
and modes that have a significant amount of their energy transported 
in the cladding are eliminated.*® Our understanding of these factors 
has largely come from approximate solutions to Maxwell’s equations 
for lightguides. The availability of computer software for exact numer- 
ical solutions’ to Maxwell’s equations allows us to reexamine some of 
these effects and to look at other factors that may influence lightguide 
bandwidths. In this paper, we determine optimum power law profile 
parameters (i.e., a parameters) as a function of wavelength for a 
germania-doped silica lightguide. 


ll. PROCEDURE 


We study a lightguide with a 13.5 mole percent germania-doped 
silica core and a pure silica cladding. The core is assumed to be 
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perfectly cylindrical with a diameter of 50 microns. The cladding is 
considered to be essentially infinite. The doping profile is such that 
the refractive index is a power law function of distance R. Thus, 





N= N+ (Na~ N-(Z), (1) 


where JN, is the refractive index at the center of the core, Ne is the 
refractive index in the cladding, R., is the radius of the core, and a is 
the profile parameter. We assume that this functional form holds at all 
wavelengths of interest. 

We have accurately determined the wavelength variation of the 
refractive index of silica and some doped silicas. Figure 1 gives the 
data for 13.5 mole percent germania-doped silica and pure silica. As is 
well known, the refractive index falls off slowly with wavelength for 
both of these glasses.® We express the data by the following modified 
Sellmeier formula:? 


N; = Co + CX + CoA‘ + C3/(? — 0.035) 
+ C,/(d? — 0.035)? + Cs/(A” — 0.035)*, (2) 
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Fig. 1—The wavelength dependence of 13.5 mole percent germania-doped silica and 
pure silica. 


456 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1981 


where A is in microns and 7 = I, 2. Table I lists the coefficients for eq. 
(2) for both 13.5 mole percent germania-doped silica and silica. 

Our numerical procedure to solve Maxwell’s equations has been 
described in detail in an earlier paper.’ Thus, we only give a very brief 
description here. The method is similar to that described by Vigants 
and Schlesinger,’ and Vassel.!! We write Maxwell’s equations in 
cylindrical coordinates as 


dT 7 
ae "A(p)I'(p), (3) 
p 
where I(p) is a column vector whose four elements are related to the 
tangential components of the electromagnetic field. They are as fol- 
lows: 


Tl, = E2, 

T, = -ipH,Zo, s 
Ds = pK, 

I, = -—tH.Zo, 


where Zp is the wave impedance of free space and the variable p is 
defined as KoR. 
The 4 X 4 matrix A(p) can be written as 


0 (N2/«) — 1 0 —MN./« 
p«—M? 0 MN. 0 

0 MN.-/« 0 o” — (M?/k) 6) 
—MN. 0 Ni-« 0 


This matrix describes the properties of the media. 

We solve eq. (3) by an optimized fourth-order Runge-Kutta proce- 
dure.” This is a higher-accuracy formula than used in our earlier work. 
The most important computational results are the effective indices N. 
for the various modes in the lightguide. 

The group indices N, are calculated from the effective indices by 
the following formula: 


Nz = N. — AdN./dx. (6) 
Table |—Coefficients for eq. (2) 


13.5 Mole Percent Germania- 


Doped Core Silica Cladding 
Co=  1.4706868 1.4508554 
C,; = —0.0026870 —0.0031268 
Cz = —0.0000356 —0.0000381 
C3;= 0.0035756 0.0030270 
C4 = —0.0000828 —0.0000779 
Cs = 0.0000018 9.0000018 
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The derivative is obtained from numerical calculations of N. at three 
closely spaced wavelengths. . 

For a lightguide with the materials and profile described earlier, a 
simple approximate relationship exists between group index N, and 
effective index N, for certain modes. In particular, for a given profile 
parameter a, angular mode number M, wavelength \, and for modes 
far from cutoff, the effective indices are nearly linearly related to the 
group indices. Evidence for this can be obtained from Fig. 2, curve (a). 
In this case, we have \ = 0.6328 micron, M = 0, and a = 2. The solid 
line is the linear least-squares fit to the data when the last two modes 
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Fig. 2—Plot of group index (Nz) versus effective index (N.) at a wavelength (A) of 
0.6328 micron with M = 0. In (a) the profile parameter a is 2.0 while in (b) a is 2.2. Note 
that a good linear relationship exists between Nz, and N- if the last few modes are 
neglected. 
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Fig. 3—Plot of group index versus effective index for A = 0.82 micron and M = 0. 
When a = 2.08 the group indices are all nearly the same except for those modes being 
influenced by the cladding. 


are deleted. (The last mode is too far off scale to be plotted.) Obviously 
a nearly linear relationship exists between N, and Ne. In addition, 
there is a very slight systematic deviation from the least-squares line. 
This suggests that additional terms are needed to completely describe 
the data. However, this will be of no consequence to the analysis to be 
described here. 

The last two modes (those that deviate strongly from the line) have 
a significant fraction of their energy being transported in the cladding. 
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Thus, they are not tightly bound to the core and normally are not 
considered in an optimum a calculation. 

An optimum profile parameter a, for a particular family of modes, 
is one that minimizes the spread in group indices for that family at a 
particular wavelength. For example, we might consider all the modes 
associated with a specified M value. In this case, the optimum a would 
be the one that gives the smallest slope to the least-squares line 
relating N. and Nz. 
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Fig. 4—Part of the data of Fig. 3 on an expanded scale. 
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Fig. 5—Plot of group index versus effective index for A = 1.06 microns and M = 0. 


Figure 2, curve (b) shows a near optimum value of a = 2.2 for M = 
0 and A = 0.6328 micron. We see that the least-squares line is now 
nearly horizontal and that evidence exists for the necessity of higher- 
order terms to exactly fit the data. The strong drop-off of the last few 
modes is still evident. 

A very good overall optimum a for all M values will be nearly 
identical to the optimum a for M = 0. As will be shown shortly, this 
results in excellent equalization of the lower-order modes and a good 
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equalization of the higher-order modes. For most engineering appli- 
cations, optimum a’s in the wavelength range of 0.63 to 1.55 microns 
are needed. 


lll. RESULTS 


Figures 3 through 7 show plots of N, versus N, for various a 
parameters and wavelengths for M = 0. Figure 4 shows part of the 
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Fig. 6-—Plot of group index versus effective index for \ = 1.32 microns and M = 0. 
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Fig. 7—Plot of group index versus effective index for \ = 1.55 microns and M = 0. 


data of Fig. 3 on an expanded scale. Note that in all a nearly linear 
relationship exists between N, and N-., except for the last couple of 
modes. As mentioned before, these modes are nearing cutoff and are 
usually not considered in an optimum a determination. 

At each wavelength, we see that there is an a value, aop:, that makes 
the slope of the line nearly zero. Table II lists a.» values for wave- 
lengths of engineering importance and also for a few other wavelengths 
as well. 

The data of Figs. 3-7 pertaining to the infrared wavelengths are 
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Table 1I—Optimum alpha values 


A (microns) Qopt 
0.6328 2.211 
0.73 2.130 
0.82 2.081 
0.90 2.037 
0.983 2.000 
1.06 1.974 
1.20 1.925 
1.32 1.884 
1.40 1.856 
1.55 1.807 


summarized in Figs. 8 and 9. Figure 8 shows the slope of the line 
relating N, to N, as a function of wavelength for various a values. We 
see that most infrared wavelengths of interest require an a value in 
the range from 1.8 to 2.1. Figure 9 shows the same data, but this time 
as a function of a for various wavelengths. The same conclusion is 
obvious. 


Figure 10 displays a calculation of the group indices for all the bound 
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Fig. 8—Plot of the slope of N, versus N, as a function of A for various «@ values. 


464 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1981 


modes in the fiber for A = 0.82 micron and aopt = 2.081. This plot 
gives a vivid representation of the fiber characteristics. We see that 
indeed a = 2.081 gives a good equalization for all the modes if those 
being influenced by the cladding are ignored. A delay-time calibration 
corresponding to 1 ns/km is also included in the plot. 


IV. DISCUSSION 


Within the framework of the WKB approximation, Marcatili’® has 
obtained an exact result for ao» that minimizes pulse widths. His 
formula is as follows: 


Qopt = (2 — p)(1 + V1 — 2A) — 2, (7) 
where 
N, A aA 
= Ne A’ am’ (8) 
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Fig. 9—Plot of the slope of N, versus N, as a function of a for various A values. 
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Fig. 10—-Plot of group index versus effective index for all the bound modes in the 
lightguide. We see that for \ = 0.82 micron and a = 2.081 there is a good equalization of 
all the group indices. 


N; — Ne 


Nis the group index on axis, and A= 7 
1 


(9) 


Also, within the same WKB framework, Olshansky and Keck” have 
obtained an approximate formula for aop: that minimizes the rms width 
of a pulse. Approximations that coincide with Olshansky and Keck’s 
can also be obtained from the exact result of Marcatili. 

It is clear that aop, depends upon the type and number of approxi- 
mations employed, how the various modes are weighted, and whether 
we wish to minimize pulse width, rms width, or some other parameter. 
Another possibility would be to define an optimum a that would 
include the effects of those modes having significant energy trans- 
ported in the cladding. This is likely to lead to an a quite different 
from that obtained in this work. 
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Figure 11 shows the results of Olshansky, Marcatili, and ourselves. 
There is great similarity in the shape of the curves and in the actual 
numbers. The deviation between our results and the WKB results is | 
largest at the longer wavelengths. It is obvious that the agreement 
between all results is good and, at the present state of the manufac- 
turing art, can be considered the same. However, this may not neces- 
sarily be the case in the future. To illustrate (see Fig. 12), we calculate 
N, as a function of N, at 1.55 microns for M = 0 and ao = 1.775 
(Olshansky) and adopt = 1.807 (this work). There is an obvious difference 
in slope; this difference would be significant in the design of high 
bandwidth lightguides. 

A good representation of our optimum a values is given by the 
following formulas: 


0.63 micron S A S 1.00 micron, 


Qopt = 2.9970857 — 1.6647237A + 0.6629181)", (10) 
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Fig. 11—Plots of optimum profile parameter a versus wavelength for a 13.5 mole 
percent germania-doped silica lightguide. The index dispersion data employed is from 


ee work. A clear similarity exists among the results of Marcatili, Olshansky, and 
this work. 
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Fig. 12—A plot of N, versus N, for \ = 1.55 microns. aopt = 1.775 is calculated from 
the formula of Olshansky et al.,.and aopt = 1.807 is from this work. 


1.00 micron < A S 1.55 microns, 
Qopt = 2.38835734 — 0.4170708A + 0.0290252A7. (11) 


It is reasonable to assume that by linear interpolation or extrapola- 
tion we can generate the refractive index dispersion curves for ger- 
mania-doped silica in the range from 0 to 16 mole percent. Using this 
data, we can then calculate the dependence of N, on N, for a number 
of levels of doping. Figure 13 shows a family of curves for M = 0, a = 
1.884, and A = 1.32 microns. There is a good equalization of N, for all 
the compositions. We also see that, as expected, the number of modes 
decreases as the doping level drops. Thus, at the 1 percent doping level 
for M = 0 only one WN, exists. As usual, the last couple of modes for 
each doping level show a substantial cladding effect. 

Finally, the data displayed in Figs. 8 and 9 may have some implica- 
tions to lightguide engineering. From Fig. 8, we see that the various a 
curves pass through zero with about the same slope. This suggests that 
the range of operating wavelengths will be the same for all design 
wavelengths. Likewise from Fig. 9, we see that the various 4 curves 
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Fig. 13—A family of plots of N, versus N, for various doping percentages of germania. 
The a value is 1.884 and M = 0. 


pass through zero with similar slopes. This suggests that the allowable 
error in a will be about the same for all design a’s. 
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Multistage Decoding 
of Frequency-Hopped FSK System 


By UZI TIMOR 
(Manuscript received July 17, 1980) 


A recent paper described an improved decoding scheme for a 
frequency-hopped multilevel Fsk system. We examined this multiple 
access communication system for possible application in satellite 
communication and mobile radio telephony. The new decoder, using 
the known algebraic structure of the users’ addresses, reduces mutual 
interference and achieves a 50 to 60 percent increase in efficiency 
over conventional decoding. The present paper shows how additional 
decoding can further increase the efficiency, bringing it very close 
(within half a percent) to optimum. The scheme makes use of infor- 
mation derived while decoding the messages of other users and thus 
is especially attractive for the base station, where such information 
is readily available and does not require a significant increase in 
complexity. Compared to conventional decoding, the new scheme 
more than doubles the number of simultaneous users. 


1. INTRODUCTION 


In a recent paper’ a new decoding scheme for a frequency-hopped 
multilevel FsK system was described. This system, where M users 
share a common frequency band, has been examined for applications 
in multiple access satellite communication” and mobile radio tele- 
phony.’ Each user conveys a K-bit message every T seconds by 
transmitting a sequence of L tones (chips) chosen from an alphabet of 
2* sinewaves of duration r(—=7/L). Each user is assigned an address 
(code) and the message is modulated onto the address. The receiver, 
knowing the address, decodes the received signal and extracts the 
message. However, transmissions by other users can combine to cause 
an erroneous message resulting in an ambiguous reception. Thus, even 
without channel impairments the number of simultaneous users the 
system can support at a given error probability is interference limited. 

Whenever an ambiguous reception occurs, i.e., more than one mes- 
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sage value is decoded, the conventional decoder can’t identify the 
correct message. The new decoder’ makes use of the algebraic structure 
of the addresses to eliminate erroneous messages that come from 
interference. Each one of the L chips comprising a message is checked 
to determine whether it is part of an interference pattern, ie., a 
possible sequence transmitted by some other user. An interference 
message will always be identified as such. The correct message will 
usually fail to have interference patterns for some chips and thus will 
be identified and correctly decoded. | 

This additional decoding (stage 2) results in a substantial improve- 
ment in performance over conventional decoding, allowing a 50 to 60 
percent increase in the number of users that can simultaneously share 
the system at a given error probability. Yet it is not optimal and it can 
be further improved. 

In performing stage 2 decoding interference patterns are sought but 
the (addresses of the) users who might have caused them do not have 
to be identified. This information is readily available and can be used 
to eliminate pseudo interference patterns. For example, an interference 
pattern might “belong” to a user who is not active at the moment. 
Even if it belongs to an active user it might not be his sequence but a 
combination of others. 

This paper describes how additional decoding (stage 3) eliminates 
pseudo interference patterns and thus further improves the perform- 
ance of the system. There is a clear distinction between decoding at 
the base station and at the mobile unit. The base station decodes the 
messages of all users and thus any information required for stage 3 
decoding is already available. Therefore, additional decoding with 
improved performance can be accomplished without a significant 
increase in complexity. In the mobile unit, on the other hand, the 
additional decoding requires an increase in complexity as the messages 
of other users have to be decoded. 

Upper bounds on the number of simultaneous users that the system 
can accommodate at a given error probability as well as simulation 
results are presented. A comparison with a hypothetical decoder, 
which fails only when the correct message coincides with an interfer- 
ence message (and therefore is undecodable), shows that the stage 3 
decoder is very close (within half a percent) to the optimum decoder. 

In the noiseless case where the system has a total bandwidth of 20 
MHz and the data rate of each user is 32 kilobit/s, the number of users 
that can simultaneously share the system at a bit error probability of 
10~° is increased from 216 (in conventional decoding) through 345 (in 
stage 2 decoding) to 450 (in stage 3 decoding)—a total improvement of 
108 percent. The resulting efficiency (total rate transmitted through 
the system per unit bandwidth) of stage 3 decoding is 72 percent. 
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Under noisy and multipath conditions, the number of simultaneous 
users is reduced but a comparable advantage of the new decoder over 
the conventional one is maintained. 


li. SYSTEM DESCRIPTION 

We briefly describe the system here. A detailed description is given 
in Refs. 2 and 3. 

The elementary signals of the system are a set of 2% sinewaves, 
which are orthogonal over the chip duration tr. Each user conveys a K- 
bit message every T' seconds by transmitting a sequence of L tones of 
duration 7 chosen from the signal set. The sequence is determined by 
the user’s address and his K-bit message. 

Let the address of the mth user be denoted by a vector an, 


Qn = (Gmi, Am2, **°; Amz); (1) 


where each am; is a K-bit number corresponding to one of the 2* 
frequencies of the system. 
The transmitted sequence is 


Ym = am + Xm-l, (2) 


where X,, is the K-bit message, and 


A 
1G 1, 24,4), 





L 


Let us denote the 2* frequencies of the system as elements of 
GF (2*), the finite field (Galois field) of 2* elements 0, 1, ---, 2% — 1. 
Accordingly, the message X,,, the components of a, and Y,, can be 
expressed as elements of GF(2*), and the operations of addition 
(subtraction) and multiplication (division) are performed according to 
the rules of GF(2*). The transmitted sequence Y,, can be described as 
a pattern in the L x 2* time-frequency matrix A. Simultaneous 
transmissions by M users will result in up to L X M entries in A. The 
receiver performs every 7 seconds a spectral analysis of the composite 
received signal and decides which of the 2* frequency cells contain 
energy. Thus, after T seconds, assuming no channel impairments, a 
duplicate of A is generated at the receiver. To decode X,,, the receiver 
adds —a,n: to each column of A to obtain user m’s decoded matrix A». 
The message appears as a complete row in A,,. Ambiguous decoding 
occurs when transmissions by other users combine to form other 
complete rows in An. 


li, ADDRESS ASSIGNMENT 


The interference between users can be minimized by a proper choice 
of addresses with minimum cross correlation. Schemes for assigning 
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2* addresses that guarantee minimum mutual interference between 2* 
or fewer users have been proposed‘ using an algebraic approach. 

Let us consider the following address structure: The address of user 
m is defined to be‘ 


an = (Ym, mB, iseetaes | mB"), (3) 


where ym is the element in GF(2*) assigned uniquely to user m and 8 
is a primitive element in GF(2*), which is fixed for the system. 

Note that the modulation (2) and address assignment (3) are not 
unique. Other schemes are possible (see for example Refs. 4 and 5) 
and the analysis and results that follow can be extended to such 
systems. 


IV. STAGE 2 DECODING 


Let M be the number of simultaneous users of the system. Each 
user 7 transmits a sequence Y; (2) where the address a; is assigned 
according to (3) and the message value is X;. 

In the decoded matrix A,, of user m, Xm appears as a complete row 
Xm:1. Suppose we have another complete row X’ in A,, that is the 
result of interference. According to the address assignment, each chip 
in X’ must come from a different user, so at least L users i, -++, iz 
have combined to cause this interference. We denote by 7, the user 
that contributes the interference in column n of X’. (If more than one 
user caused this interference, i, can be any one of them.) To simplify 
notation, let us denote the address element (y;,) and the message value 
(X;,) of user in by yin and Xin, respectively. 

As mentioned in the previous section, the number (location) of each 
row in A,, is an element of GF(2"). We can subtract (the number of) 
row X’ from all rows in A,, to obtain a new matrix Dx:, where row X’ 
will now be at row zero. 

Let gn(7) denote an entry in row q, and column j in Dx:. 

It was shown’ that a necessary condition for row X’ in A, to be 
caused by interference is the existence of L nonzero numbers 6, in 
GF (2*) such that all entries, 


Qn(J) = 5n(B7? Br), ds n= 1; ey L, (4) 


appear in Dx:. The entries g,(j) are the contribution of user i, with 
address element yin = 5, + Ym, where ym is the address element of user 
m. From (4) we can derive the following relation: 


Qn(n + 1) = Qnl(J) fi-n; n=1,---,L-1, (5) 
qu(L — 1) = qi(J) fF2, (6) 
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where 


fn pom, Fe L, Jen, (7) 
ff-.= = fi-L. (8) 
B 
To check the condition (4) for n (n = 1, ---, L) we multiply all entries 
(i.e., their row numbers) of column] in Dy: (j= 1,---,L,j#njA#An 


+1) by f-n (or fz) and look for a “complete” row (complete except 
the nth term). If we find such a row, then there exists a possible 
interferer contributing to column n of X’. 

If the condition is satisfied for all n = 1, ---, L, we assume that row 
X’in A, is the result of interference. Using this decoding scheme, all 
interference rows will be identified as such. The correct row X,, will 
usually fail to satisfy the condition for some n and thus could be 
identified and decoded as the correct message. 

It was shown’ that a simple way to perform the column multiplica- 
tions required during the test is to express each row number by an 
exponent of 8. Since 8 is a primitive element of GF(2*) this is 
equivalent to a (fixed) row permutation of Dx. This way the multi- 
plications can be substituted by cycle shifting of columns, which is 
easy to implement. Thus, stage 2 decoding requires a small increase in 
complexity (of an order of L’ cycle shifts) yet achieves a substantial 
improvement in performance. We show below how additional decoding 
can further improve the performance. 


V. STAGE 3 DECODING 
5.1 Discussion 


Stage 2 decoding fails to identify the correct message X,, when for 
every n (n = 1, ---, L) a possible interferer 1, exists such that all 
entries gn(j),j = 1, ---, L, appear in Dx. 

If such L interferers actually exist, that is, if X,, coincides with an 
interference row, there is no way to distinguish between X,, and any 
other interference row and no decoder can identify the correct message. 

There are however many cases when only L — / chips of X,, are 
(also) the result of interference and the remaining / chips, although 
not caused by interference, satisfy the interference condition. That is, 
entries from other users combine to form interference patterns at those 
chips. If those pseudo interferences can be identified, for at least one 
chip, X, will be identified and correctly decoded. This will be done by 
making use of available information, derived but not used in stage 2 
decoding. 

Recall that a possible interference pattern at chip n of X,, was 
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identified as a “complete” row (complete except for the term in the 
nth column) in Dx. From (4) to (6) this row will be at row number q,, 
where 


_ Jd(B"-B), n=1,-+-,L-1, 
Qn = tae. = Be), n=L. (9) 


Knowing gq, (which is found during the test), we can easily compute 
6, and the address yin = 5, + ym of the possible interferer. The message 
value X* associated with this interference can be computed as follows: 
The transmitted sequence of user 7, with message value X™* is (2), 


Yin, = Yin? + X*.1. (10) 
The entry at column n of A,, will be at row number 
YinB”* + X* — ymB"~* = bnB" + X*. (11) 


But by definition i, is the user who contributes the interference at 
column n of X,,. Thus, 


5,8") + X* = Xm 


or 
Ke SX O, 0" (12) 
Substituting 6, from (9), we can write 
X,- i ee Pe 
X*t= a (13) 
Xm —-———, n=L. 
L= £8 


5.2 Principle of decoding 


Let us assume that the list of (the addresses of) the active users 
(those who are currently using the system) is known. This is certainly 
the case for the base station that communicates with all its active 
users, but the list could also be relayed to the mobile unit. 

Thus, if an address y;, of a possible interferer belongs to a user who 
is not active, the interference can be immediately identified as pseudo 
interference and eliminated. 

If the address y;, belongs to an active user, there are two possibilities: 

(1) X* = Xin. The transmitted message Xi, and the computed 
message X * associated with the interference pattern are the same. In 
this case the interference was actually caused by user ip. 

(tz) X* # Xin. In this case the interference pattern was not caused 
by user z,, but is the result of a combination of (at least) Z — 1 other 
users. Thus, this is a pseudo interference and should be eliminated. 

Consider the decoded matrix A;, of user 7,. His transmitted sequence 
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will appear as a complete row at row number Xn. If X* ~ Xin, then we 
will have another complete row (an interference row) in Aj, at row 
number X*. By decoding Ain (stage 2 or stage 3 when necessary) the 
correct message Xin will usually be identified. Thus, if X* 4 Xin, it will 
be identified as pseudo interference of user m and discarded. 

If there are several possible interferers at column n of X,, (more than 
one “complete” row is found in the test) all of them must be checked 
and discarded before we can conclude that interference is not at that 
column. If we find (at least) one of the L chips for which no actual 
interference exists, we can conclude that X,, is the correct message. 

If we cannot assume knowledge of the list of active users, then the 
pseudo interference pattern attributed to nonactive users cannot be 
identified as such since there is no message X;, to be decoded. 


VI. DECODING AT THE BASE STATION AND AT THE MOBILE UNIT 
6.1 Base station 


The base station has a list of all active users and decodes their 
messages. The decoding procedure can be described as follows. 


Step 1 (for each user m) 


(t) The decoded matrix A,, is generated and searched for complete 
rows. If there is only one such row, the message is decoded and listed 
in a list (L1) of users whose decoding is completed. 

(tt) If there is more than one complete row, stage 2 decoding is 
performed. If the message is decoded it is transferred to L1. 

(tit) If the message is not identified, the address and message values 
of possible interferers are computed [according to (9) and (13)]. Those 
are compared with the list of nonactive users and with L1 to eliminate 
a pseudo interference pattern. 

At any stage, if the message is decoded it is transferred to L1. 
Similarly, if the message is found to be undecodable (i.e., all interfer- 
ence patterns are found to be true interferences) it is transferred to list 
L2 of undecodable users. 

At the end of step 1 (for all users) we have N; users in list L1, No 
users in L2 and N3 = M — (N, + Ne) users, each with a list of possible 
interferers and associated message values for some of its chips. If Ns 
> 0, we proceed to step 2. 


Step 2 (for users who have not been decoded) 


For each of the N3 users, we compare his list of possible interferers 
with L1 to eliminate pseudo interferers. As list L1 is increased, more 
possible interferers can be checked and more users decoded. 

The procedure terminates when N3 = 0 or when L1 remains un- 
changed after a complete cycle of checking all Ns users. 
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6.2 Mobile unit (user m) 


(t) The decoded matrix A,, is generated and searched for complete 
rows. If (up to) stage 2 decoding does not identify the message a list of 
possible interferers and their associated message values is computed. 

(it) Pseudo interferences from nonactive users (if list is available) 
are identified and discarded. 

(uit) If the message is still undecoded, the decoded matrices of 
possible interferers have to be generated and decoded (stage 2) until 
enough pseudo interferences are identified to decode the message. 

We define this as stage 3/2 decoding since the messages of possible 
interferers are decoded up to stage 2. 

We could in principle perform stage 3 decoding of the possible 
interferers; however, the complexity increases exponentially as more 
and more users have to be decoded. 


7. ERROR PROBABILITIES 
7.1 Upper bounds 


We have an ambiguity in decoding when the following two conditions 
are satisfied. 

(t) There are two or more complete rows in An. 

(ti) The correct row cannot be identified, i.e., all its chips have 
interference patterns (actual interference or unidentified pseudo inter- 
ference). 

When this happens, we choose one of the complete rows at random 
and decode it as the message. 

Although not strictly independent, it can be shown using random 
coding arguments that the two conditions can be assumed to be 
independent with a negligible effect on the probability of error. Thus, 
PP, the word-error probability when stage i decoding is performed, is 
given by 


P® = P,P#, (14) 


where P; is the probability of condition 1 and P$? is that of condition 
2 when stage 7 decoding is performed. 

If 2* is the number of frequencies, L is the length of the sequence 
and M is the number of simultaneous users, P: can be upper-bounded? 
by 


P, < (2% — 1) p’, (15) 

where 
pHleadag yes (16) 
Let P§; be the probability of condition 2 (stage i decoding) when j 
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chips in X,, come from interference and the remaining L-j chips have 
a pseudo interference pattern. We then have 


L 
PS = Y P§,. (17) 


j=0 


(a) Stage 2 decoding 


For stage 2 decoding we have shown’ that P’ can be upper-bounded 
by 


Pr< ({ Js (18) 

where 
S = (2* —1)(1—p)p*? (19) 

Thus, 
P? <p*(14+S)*, (20) 


and as long as p(1 + S) < 1 the word-error probability can be upper- 
bounded by 


P® < (2% — 1)(1 + S)*p” (21) 
and the bit-error probability by 
PP < 2*-7(1 + S)"p™" (22) 


(b) Stage 3 decoding 


P$), corresponds to the case where all interference patterns actually 
come from other users (i.e., X,, coincides with an interference row in 
Am). Thus, 


Pot = Py <p" (23) 


Consider the case where L-1 chips in X,, actually come from interfer- 
ence and the remaining chip has a pseudo interference pattern. This 
pseudo interference will correspond to an active user with probability. 


(24) 
The probability that the correct message of this user (and therefore 


the pseudo interference) will not be identified as such when stage 3 
decoding is performed for that user is P§’. Thus we can upper-bound 


Poa < Po PP. (25) 
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Similarly, 


P23 < PEA PPP’. (26) 
From (17), (18), and (26) we get 
P¥) < p*(1 + SpP¥)-. (27) 


Let P# be the solution of (27) when an equality is substituted for the 
inequality. Then 


PP < PE = p"(1 + SpP3)". (28) 


(c) Stage 3/2 decoding 


In a similar way we can upper-bound P$”, the probability of error 
of condition 2 when the messages of possible interfers are decoded up 
to stage 2. This is relevant only to the mobile unit (see section 6.2). 


Case 1: List of active users available 


Substituting P?) for P® in (25) yields 


PY” < p*(1 + SpP?)’, (29) 
where P is given by (20). 


Case 2: List of active users unknown 


In this case we will always have an error when a possible interferer 
is not an active user (which occurs with probability 1 — p) since there 
is no correct message to be decoded. Thus 


PS” < p*[1+ S(pP? +1-p)]” (30) 


(d) Coincidence of X,, with an interference 


Whenever we have two or more interference rows in A, and the 
correct row coincides with one of them, it cannot be identified and 
decoded. An optimal decoder can do no better than a decoder which 
fails only when such a coincidence occurs. Thus, it is of interest to 
compare the various schemes with such a (hypothetical, not necessarily 
achievable) decoder. 

The probability of condition 2 for such a decoder is given by (23) 


PE < p*. (31) 


We can summarize the performance of the various schemes as follows: 
The word-error and the bit-error probabilities are upper-bounded by 


PY < (2% — 1)gip”, (32) 
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PO arp, (33) 


where 
1+ 5S, stage 2 decoding, 
1+ SpPi, stage 3 decoding, 
gi= 11+ SpP, stage 3/2 decoding 
(list of users known), (34) 


1+ S(pPY+1—p), stage 3/2 decoding 
(list of users unknown), 
1, coincidence fail decoder, 


and S, P¥, p, and P# are given by (19), (20), (22), and (28), respectively. 

The upper bounds of the various schemes as a function of the 
number of simultaneous users (M) are depicted in Fig. 1 for the case 
K = 9 (512 frequencies) and L = 11 (which is optimal when the total 
bandwidth is 20 MHz and each user’s rate is 32 kilobit/s). 

As can be seen, at small error probabilities the performance of stage 
3 decoding is very close to the probability of coincidence, that is, 
almost all messages which do not coincide with an interference are 
correctly decoded. As the number of users is increased, a strong 
threshold effect occurs (at M = 475 in the above example) and the 
performance deteriorates rapidly. This is due to the fact that as M 
increases, the interference is so high that a large number of users 
cannot be decoded (stage 2) and thus the information (decoded mes- 
sage values of possible interferers) required for stage 3 decoding is not 
available. 

As M is decreased, all the above bounds converge to the coincidence 
bound. 


7.2 Simulations 


The various decoding schemes were simulated on a Digital Equip- 
ment Corporation PDP 11 computer. Figure 1 shows simulation results 
for stage 2, stage 3, and stage 3/2 (when the list of active users is 
known) decoding as well as for the probability of coincidence. 

It can be seen that for bit-error probabilities up to 107%, the simulated 
results are within 5 percent of the upper bounds on the number of 
simultaneous users. 

The threshold effect of stage 3 decoding is clearly seen. As M 
increases there are some K-bit blocks where the majority of users (60 
percent and more) are undecodable. The percentage of undecodable 
blocks increases from 6 percent at M = 478 to over 75 percent at M 
= 490. 
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Fig. 1—Upper bounds and simulations of the bit-error probabilities as a function of 
the number of users M for stage 2, stage 3/2 (list of active user known), and stage 3 
decoding and for coincidence of the message with interference. The system parameters 
are W = 20 MHz, R = 82 kbit/s, K = 9, and L = 11. 


Vill. CONCLUSIONS 


The performance of the frequency-hopped multilevel Fsk system 
can be substantially improved by making use of the known structure 
of the addresses to perform additional decoding. In stage 2 decoding 
additional checks are performed on the user’s own decoded matrix but 
no reference to the decoded matrices of other users is needed. The 
stage 3 decoder achieves a performance that is very close to optimum 
by referring to results obtained in decoding the sequences of other 
users. 
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In the base station, where all users are decoded anyway, stage 3 
decoding can be realized without a significant increase in complexity. 
This is not the case at the mobile unit; therefore, an intermediate stage 
of decoding (stage 3/2) was proposed, requiring a moderate increase in 
complexity. 

When the total bandwidth is 20 MHz and the transmission rate of 
each user is 32 kilobit/s, the number of simultaneous users that the 
system can accommodate at bit-error probability of 10~° is increased 
from 216 (conventional decoding) through 345 (stage 2) to 383 (stage 
3/2) and 450 (stage 3). The corresponding efficiencies of the system 
(total rate transmitted through the system per unit bandwidth) are 35 
percent (conventional), 55 percent (stage 2), 61 percent (stage 3/2), 
and 72 percent (stage 3). As can be seen, stage 3 decoding more than 
doubles the efficiency. Simulation results show that an efficiency up to 
75 percent can be obtained with stage 3 decoding. 
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A Bubble Memory Differential Detector 


By W. D. WYNN 
(Manuscript received September 5, 1980) 


In the application of digital memories, significant variations can 
occur in the memory sensor analog responses processed by detection 
circuitry to recover stored binary data. Across an ensemble of memory 
devices of the same type, the variations are due to manufacturing 
deviations, and for a particular device, they occur with changes in 
operating conditions. With magnetic bubble memories, binary data 
are recovered by sensing the presence or absence of a magnetic 
domain in each bit interval. In general, a nonadaptive bubble detector 
does not accommodate the variations that may occur in the bubble 
response when each decision is based only on the sensor response in 
the current bit interval. Here we present a differential detector design 
that is reasonably simple to implement as one integrated circuit and 
is tolerant of signal amplitude drifts in bubble memory sensor outputs. 


1. INTRODUCTION 


In typical field-accessed magnetic bubble memories (MBM), stored 
binary data on memory chips are detected by processing millivolt-level 
signals that are measured differentially across a matched pair of 
permalloy magnetoresistive sensor strips.’* Ideally these sensors are 
subjected to the same magnetic, electrical, and environmental condi- 
tions with the exception that only one, the active sensor, is subjected 
to H fields of traversing bubble domains. In the drive field period of a 
field-accessed MBM chip, by design at least one time interval exists, 
defined here as the detection interval (D1), in which the sensor differ- 
ential response is relatively free of unavoidable signal interferences 
such as MBM function drive crosstalk and sensor magnetoresistive 
switching noise. Within the DI a detectable difference should exist 
between responses for the bubble and no-bubble cycles of the memory. 
Even after carefully choosing the DI, undesirable sensor response 
variations can still occur. Variations within the pi that have been of 
major concern in the use of previous nonadaptive MBM detectors are 
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the device processing variation in magnetoresistive sensor response, 
the negative temperature coefficient of this response, the phase shift 
of bubble response transition that occurs with chip manufacturing and 
drive field changes, and the adjacent symbol interaction of bubbles in 
the sensor array. In a detector where the pi is fixed and where the 
individual bubble response in the DI is compared to a fixed threshold 
these variations have been intolerable. 

An accurate analytical parametric description of the family of sensor 
signals that are observed for any MBM chip type has been intractable 
because of the nonlinear field effects of permalloy magnetoresistive 
sensor arrays and a lack of models for statistical data about the 
manufacturing variations for MBM devices. Figure 1 shows qualitative 
examples of differenced sensor responses within a DI that may exist for 
a chip type when the manufacturing and operational variations are 
considered. Each signal pair is shown in the DI for some particular chip 
of an ensemble of the same chip type at a particular operating condi- 
tion. All signals have been referenced to a voltage V, at the start of the 
DI. The amplitude and time scales for Fig. 1 are representative of field- 
accessed MBMs with field drive frequencies in the 100-kHz range. Not 
only can response amplitude change with operating conditions and 
selected chip, but also the shape of the desired bubble response and 
intersymbol interference contributions may vary. 

Detectors could be designed that adapt to variations in the bubble 
and no-bubble responses where such variations preclude a fixed thresh- 
old detector. For example, detector amplifier gains could be tempera- 
ture compensated to track the temperature sensitivity of the magne- 
toresistive sensor. Also, the phase of the DI could be adaptive, and 
detection thresholds could be adjusted dynamically from responses in 
a known data pattern. Unfortunately, adaptive bubble detectors have 
disadvantages such as increased circuit complexity, increased cost, and 
system limitations such as detector adaptation time before reading 
data. 
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Fig. 1—Dc restored bubble (denoted by “1’’) and no-bubble (denoted by “0” ) differ- 
ential sensor responses for a memory chip ensemble. 
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We describe a differential detection approach, and a circuit realiza- 
tion of it, which are insensitive to some signal amplitude, phase and 
shape variations in the pI, if these variations occur slowly when 
compared with the memory drive field period. The detector is based 
on the hypothesis that some signal function, such as peak-to-peak 
transition or the signal energy, is available in a DI such that the 
minimum functional change between unlike data pairs is greater than 
the maximum functional change between like data pairs in memory 
detection cycles that are contiguous. This condition must hold over 
the range of manufacturing and environmental changes in the sensor 
outputs for the chip ensemble, including noise and intersymbol inter- 
ference effects. This proposed detector can be developed as a small 
inexpensive Ic. It is easily incorporated in MBM circuit designs, and it 
is compatible with memory system requirements. 


ll. DIFFERENTIAL MBM SIGNAL DETECTION 
2.1 Signal and interference characteristics of memory output 


As indicated by the qualitative results shown in Fig. 1, the dynamic 
range of bubble and no-bubble responses in the DI may be significant 
across the desired operating range for an ensemble of a MBM chip type. 
However, for a given MBM chip and stationary operating environment, 
the responses for two data 1’s (bubbles) or for two data 0’s in successive 
pi’s differ only through random noise and intersymbol interference 
(1s1) effects. 

Sensor responses for bubble and no-bubble traversals depend on the 
memory sensor design, the chip technology, memory manufacturing 
quality control, and operating environment. Typically, the bubble 
response transitions of the MBM sensor in the Di have been in the 2- 
mV to 20-mV range, with a more restricted dynamic range over the 
ensemble and operating range of any particular memory type. Al- 
though unit-to-unit variations may exist in both bubble and no-bubble 
responses, as indicated by Fig. 1, a typical memory product gives a 
minimum transition difference of about 2 mV between the two re- 
sponses in the DI. 

When the phase of the pI is properly selected relative to sporadic 
random interference, such as magnetoresistive switching noise from 
the sensor or control signal cross talk, residual noise within the DI 
becomes a second-order effect. For example, the important frequency 
band of bubble signals in the D1 is below 10 ff, where f; is the MBM field 
drive frequency. For MBM drive frequencies of less than 200 kHz, the 
signal processing band for a detector is then less than 2 MHz. For a 2- 
MHz noise equivalent bandwidth in a detector preamplifier, the total 
random noise within the D1 is less than 50 pV rms referred to the 
output contacts of the MBM chip sensor. This includes random noise 
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from the MBM sensors under active memory conditions, the sensor bias 
current sources, and the detector electronics. 

Fringe field coupling of bubble domains can exist across elements in 
a sensor array. This effect is observed as ISI in the DI response. Like 
the desired DI response, the IsI depends on the sensor array design 
and the memory operating conditions, but it is not well understood. 
Because of the nonlinear magnetic interaction of permalloy sensor 
arrays and bubble domains, accurate modeling and analysis of IsI as 
well as the desired bubble response have been intractable. In memory 
design, arrays have been selected to reduce the ISI, as well as random 
noise, to acceptable levels relative to the minimum difference between 
isolated bubble and no-bubble transitions in the pI. This condition is 
indicated by Fig. 2a for contiguous Dis for two signal amplitude cases. 
If Ist is severe for a memory design, this design can be changed or 
signal equalization can be considered for the detector.*” However, 
because of the variable nature of IsI as well as the desired signal’s 
shape in the DI, an equalizer should be adaptive. For typical useful DI 
durations of a few ps and very short setup times of the memory system, 
adaptive equalization of Ist is not cost effective in a bubble memory 
detector. By careful sensor design and processing control for the 
memory, it has been practical to hold Ist and noise to levels that 
permit differential detection with probability of error Pz < 10~° with- 
out using ISI equalization. 


2.2 Data precoding and differential recovery from the MBM 


Figure 1 indicates that variations in the absolute levels, phases, and 
shapes of bubble and no-bubble responses along with a close match of 
like symbol responses in contiguous Dis make differential amplitude 
detection attractive for field-accessed MBM. This is analogous to dif- 
ferential phase detection used in communication channels.° If the 
binary data “1” is stored in the MBM chip as a bubble-to-no-bubble 
transition, or the reverse transition in contiguous memory cycles, then 
recovering data from memory is possible when the absolute values of 
the greatest DI response differences for pairs of like symbols are less 
than the smallest absolute values of differences for unlike symbol 
transitions. Figure 2b summarizes these prerequisite conditions for 
successful difference detection. 

Direct difference encoding of the original data has the potential for 
propagating errors since a single error in decision can result in the loss 
of phase reference and cause ambiguity in the detector. This problem 
is eliminated by a simple precoding step for the data prior to recording. 
The precoding does not reduce the data throughput or storage capacity 
of the memory. 

An original data sequence defined by {a,}, where a; € {0, 1} for k 
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Fig. 2—(a) Dec restored bubble sensor response examples showing effects of dynamic range, IsI, and noise. (b) Response peak 
difference C; ranges and I’, [_ threshold placement. 


= 0, is to be recorded as transitions between bubble and no-bubble 
states in consecutive field cycles in a way that eliminates propagating 
detection errors through loss of reference. This is accomplished by 
precoding {a;} into an equivalent sequence {b;}, where b; € {0, 1} for 
k = 0, by the algorithm 


bp = an ® Bp-1, where 6_; = 0, (1) 


and @ denotes modulo 2 addition. Then {0;} is recorded in the memory 
with a bubble state for 6, = 1 and a no-bubble state for 6, = 0. The 
precoding in (1) is implemented with a shift register and an EX-OR 
logic stage in the memory recording circuit. 

If, in response to b;, a functional B, is available from the response 
in the DI and a decision of bubble or no-bubble in cycle k& is to be 
based on B,, difference detection is performed as follows. In the kth 
data cycle the analog difference C, 4 B, — Bz-1 is formed. Ideally, 
there are three possible analog levels for C;, corresponding to the four 
permutations of (Bz-1, B,), and a ternary decision on C;, results in a 
three-state logic output c, € {—1, 0, +1}. In terms of the binary 
sequence { b;}, 


ce = bp — bra, k=0. (2) 


From expression (1), the modulo 2 subtraction of b,-; from both sides 
of (2) gives 


(Cz)mod2 = (bz — bg-1)moa2 = Ah. (3) 


It follows then that {a;} is available by simple logic operations on the 
ternary sequence {c,} obtained from {C;,} without the problem of 
error propagation due to phase ambiguity if a detection error occurs. 

In practice, the differential detector outputs are estimates of the 
components of {c;,} that are based on the available analog differences 
{C;,}. Then the detector outputs form the estimate sequence {¢,} and 
the estimate of the original data sequence is 


{Gz} = {(Cz)moa2}, k=O. (4) 


If in the kth cycle there is no detection error, ¢;, = c, and hence a; = 
a,. Figure 3 gives an example of the operations discussed in this 
section, where the functionals {B,} are the positive peaks of dc 
restored bubble and no-bubble responses in the detection intervals. 


2.3 Choice of functional B,, 


As discussed above, many variations can occur in the sensor re- 
sponses from a given type of MBM chip which effect the shape, 
amplitude, noise content, and position of the responses in a DI. To 
successfully detect differences, any signal processing of the responses 
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Fig. 3—Differential detection example. 


leading to a sequence { B;} must result in nonoverlapping amplitude 
bands for the difference functionals { C;,}, as shown in Fig. 2b, in spite 
of the inherent variations in bubble sensor response. A diverse selection 
of practical approaches in signal processing can be considered to obtain 
{ B,} from DI responses. Examples of processing over the pI for B; are 
peak-to-peak acquisition, squaring-with-integration after dc restore, 
matched filtering of the nominal bubble response, and peak acquisition 
after dc restore. Each approach has advantages and disadvantages in 
detection error performance and detector circuit complexity for field 
accessed MBM signals as understood to date. Because of the difficulty 
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in modeling the statistical characteristics of the MBM signal and inter- 
ference components, the analytical solution for a useful optimum 
detector has been intractable. The selection of any detector, and, in 
particular any signal processing approach to obtain {B,} for the 
differential detector considered here, has been based primarily on the 
performance observed in experiments with a representative ensemble 
of devices for each MBM design. The optimum choice of a detector is 
complicated additionally by practical considerations of its circuit re- 
alization, cost, and flexibility of operation. 

Section 3 describes a differential detector that has been successful 
with signal characteristics of present bubble memories. This realization 
uses peak detection of dc restored signals in the D1 to form {B,}. We 
have found that signal equalization to reduce ISI and special noise 
filtering are not necessary in this detector for the bubble memory 
families of interest. 


2.4 Detection errors and monitoring 


In Section 2.2 the detected data sequence is 


{Gn} = {(Ck)moa2}, (5) 


where ¢; is the ternary logic sequence obtained by a bipolar threshold 
comparison of C, = B, — Bz-1. With the exception of noise and 
intersymbol effects in the generation of {B,}, only two values occur 
for B, under the conditions of static operation between memory 
detection cycles. Then of the nine combinations of cz, cz41, the pairs. 
(1, 1) and (—1, —1) can occur only through detection error. The 
sequence {¢; } can be monitored for the disallowed patterns as a simple 
error check, but this is not a complete error test.”* Detection errors 
that yield valid cy, cr+1 pairs cannot be detected by this simple 
monitoring. In general, an error within the memory that results in a 
recovered bubble sequence other than the {0,} that was recorded 
cannot be detected with this pattern monitoring. The altered bubble 
sequence in the MBM appears at the detector as a different {b,} 
sequence but one that could be generated by an admissible {a,} 
memory sequence. 

Since the difference detector considered here involves the threshold 
comparison of C;, = B,; — Bx-1, it would seem that a detection error in 
recovery of az = (Cz)moa2 caused by a disturbance in acquiring a B;, 
would be highly correlated with a detection error in @p+1 = (C41) mod 2 
since Cri: = Br+i — By. This apparent coupling exists regardless of the 
data formatting discussed in Section 2.2 that is used for preventing 
error strings caused by detector reference loss. Whether it is more 
probable in difference detection to have an error pair given that one 
error occurs will depend on the method of generating functional B; 
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and the interference that causes detection errors. It is known that in 
the case of additive Gaussian noise and matched filtering of signals 
without intersymbol interference, differential detection does not en- 
hance the probability of double errors given single errors.’ 


lit. A PARTICULAR DETECTOR CIRCUIT REALIZATION 


The circuit shown in Fig. 4 performs the functions discussed in 
Section 2.2, including forming the functional C, = B,; — Br-1 required 
in difference detection. The system characteristics outlined below are 
compatible with available monolithic bipolar 1c technology. Speed 
requirements implied by memory field frequencies up to 200 kHz seem 
to make bipolar technology essential at present for any analog differ- 
ence detector type of 1c. The realization of the functional C;, is simple; 
it requires a capacitor, a clamping switch, and a buffer amplifier. Then 
to obtain C; it is not necessary to construct a true sample-and-hold for 
By-1, which is an advantage in bipolar circuit integration. All capacitors 
can be integrated and the complete bubble detector implemented as a 
single inexpensive Ic with three synchronous input logic controls. 

Figure 5 shows the key voltage waveforms for the circuit for the 
example data sequence {a,} in Fig. 3. To simplify the operating 
description, any voltage offsets in switches SW1, SW2, and SW3 are 
neglected as are the droop effects due to leakage currents from capac- 
itors C2, Cs, and C, when these are in their hold states. Adverse effects 
of switch offsets in an Ic realization of the circuit can be eliminated by 
component trims during manufacture. At points in the circuit where 
leakage currents from capacitors are critical to performance, bipolar 
field effect transistor (BIFET) amplifier isolation is assumed; these are 
shown in Fig. 4. BIFET amplifiers have bias requirements on the order 
of 1 nA within a circuit ambient range of [—10°C, + 70°C]. 

The data cycle period in Fig. 5 is T = 1/f;, where f; is the MBM drive 
field frequency. The signal detection interval.(D1) is [t, t2], modulo T. 
The other intervals of interest are [ts, t,], modulo T, where switch 
SW2 is closed, and [#5, te], modulo T, where C3 is discharged through 
SW3. The maximum operating field frequency is limited by the bubble 
memory technology. The period T considered for the differential 
detector design is in the range [5 ps, 40 us]. 

In each pi, the B, acquired is the maximum of the peak catcher 
input 


ef =e, -— C2 (t1) a V;, for t€ [d, ts], modulo T, (6) 


where €.2(¢,) is the initial condition of the voltage across C2 at the 
start of the DI, and V, is the switch voltage reference. The value of 
€.2(t1) is a function of ei(t), for ¢ = t, as well as the filter action 
provided by the Roe, C2, Rs: circuit when SW1 is closed. The noise that 
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is present with the desired signal input e, of the dc Restore I circuit is 
enhanced in the restored signal ef. Because of this, a nonzero time 
constant, ta = (Roz + Rs) C2, can be advantageous. If ta — 0, then 
€.2(t1) > e1(t:), and ef (t;) — V,, which gives an exact restoration of 
ef to reference V, at ¢,. But then the steady state noise variance of ef 
in DI is twice the noise variance of é. 

The functional B; is retained on the storage capacitor C3 until 
discharged by a HI state of E3 during [ts, te]. With Z2 held LOW during 
the pi as well as the PEAK RESET interval, the Restore Switch II 
output eF follows the acquisition of B, and the reset transition that 
occurs in the peak catcher output e, when C3 is reset by SW3. When 
E2 is HI, SW2 is closed, and a voltage e,4 is stored across C, that is a 
function of both e2(t) for t € [ts, t,] and the filter time constant 7, = 
(Ros + Rs2) Cy. Then for ¢ € [t4, T'+ ts], modulo T, SW2 is opened, and 
the comparator (+) input is 


es = €2 — €ca(ts) + V,. (7) 


If t, — 0, then e.4(t4) — e2(t4) and an exact restoration of ef to V, 
occurs. But if noise is present in és, the filtering provided when 7; 4 0 
may result in improved performance to a degree that depends on the 
signal and noise characteristics of és. 

For the time sequence of controls, Ei, Ez, and E3 shown in Fig. 5, 
the value of ef(t) when ¢ = (RT + ¢3) and k = 0, 1, 2,---, is AsC, = 
A3(B, — Bz-1), where As is the gain of the second interstage amplifier. 
When, in response to like and unlike symbols in the data pair (b;-i, 
bz), the analog values C, = B, — B,-1 lie in three nonoverlapping 
bands, as shown in Fig. 2b, the recovery of the original sequence {az} 
follows from the trilevel comparison of each A3C; value with a thresh- 
old pair (T_, T+). If AsC, = T+ = V; + ys, the trilevel logic response is 
cr = +1; if AsC, = T_ = V; — y-, the decision is c, = —1; and for all 
other values of A3C;, cz = 0. Then from eq. (3), the original sequence 
{a,} is the output of the EX-NOR gate in Fig. 4. Amplifier A; in 
conjunction with the dc Restore II circuit provides the desired differ- 
encing C; in each detection cycle. A complete sample-and-hold func- 
tion is unnecessary since the difference C; and not the absolute levels 
of B, and B,-1 are required for differential detection. 

In the kth data cycle, a; is latched in FF1 by E2. The same positive 
edge of E, clocks FF2 to latch the logic level out of comparator CP1. 
The DATA and DATA MONITOR outputs are sufficient to detect 
invalid (cz, Cr+i1) pairs (+1, +1) and (—1, —1) that may occur if there 
are detection errors. This was described in Section 2.4. 


IV. DETECTOR OPERATION EXAMPLE 


The detector shown in Fig. 4 was constructed from discrete BIFET 
and bipolar transistor devices that are typical of those that could be 
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used in making a monolithic linear 1c. Table I lists a summary of the 
conditions for this discrete component circuit. We chose the param- 
eters listed in Table I, including values of the two thresholds I, and 
I_, for an MBM type that seems to be typical of the field-accessed 
bubble-memory art. For other MBM device designs, it may be desirable 
to have different parameter values than those considered here to 
maximize detection performance. 

Fig. 6a shows superimposed bubble and no-bubble responses for e; 
and the restored signals ef and ef for a recorded {b,} sequence in the 
MBM that is the periodic hexadecimal pattern ABEF5410. One of the 
two orthogonal coil currents, J,, that comprise the bubble memory 
square wave drive field, is shown also. In the ef trace of Fig. 6a, the 
detector “EYE” diagram is formed by the three levels at the right side 
of the peaks of ej. The levels of this EYE correspond to the trilevel 
logic outputs c, € {—1, 0, +1}. The thresholds [, and IT, at the 
negative inputs of comparators CP1 and CP2, are positioned within 
the EYE to get minimum detection error probability over the memory 
ensemble and its range of operating conditions. Nonadaptive detection 
is possible if '.. and I can be placed between the center band and the 
upper and lower bands of the ef EYE. The difference detector per- 
formance is immune to the absolute levels of the peaks ofey that occur 
within the detection intervals, but it does depend on the minimum 
change in peaks for unlike contiguous data pairs and the maximum 
change for like data pairs. The widths of the EYE are determined by 


Table |—Detector Parameters 


. Amplifier As: gain = 3, frequency response same as Ao. 
. Capacitors: Ci, C2, Cs, Cy = 50 pF. 
. Resistors: 
Ro, Rop = 200 Q, 
Roo, Ros = 10° Q, 
Ra, Ra, Rss = 200 Q, 
R, = 10° Q. 
10. Thresholds*: 
Tl, = V,+ y. + [closed SW2 offset], 
T_ = V,-— y- + [closed SW2 offset], 
y+, y- = 345 mV. 
11. Voltages: V, = 10V, V, = 4V, V.c = 12V 
12. Leakage currents: For ef,e,, and e3 nodes, current leakages are approximately 1 nA. 


1. Sensor resistances: Rp, Rpp = 1600 Q 

2. Bias currents: Im, Ig2 = 3.75 mA 

3. SW1 filter time constant: [ Roz + Rsi}]C2 = 60 ns 

4. SW2 filter time constant: [ Ros + Rs2|C, = 60 ns 

5. Amplifier Aj: gain = 10, 3 dB frequency = 3 MHz, critically-damped two-pole 
response. 

6. Amplifier Ay: gain = 5, 3 dB frequency = 3 MHz, two-pole response with 1 dB 
overshoot. 

7 

8 

9 


* The thresholds , and T_ are offsets from the reference voltage { V, + [Closed 
SW2 Offset}} of the comparator input ef. This is shown in Figs. 4, 5, and 6. The optimum 
values for the offsets, y+ and y-, are determined by the particular properties of the 
bubble memory family for which differential detection is considered. For the memory 
code under investigation here, the values selected were y+ = y- = 345 mV. 
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the difference between bubble and no-bubble signal transitions in the 
DI and by intersymbol interference and noise. The effect of noise has 
been of secondary importance when DI is properly selected. Also, 
through proper MBM sensor design, ISI effects due to bubble coupling 
in the sensor array can be held to an acceptable fraction of the 
minimum isolated bubble and no-bubble difference. The signal, noise, 
and IsI are all functions of the MBM design. 

Figure 6b shows the timing sequence of FE, E2, and E3 relative to 
the response ef out of the Restore Switch II for the sequence {0,} at 
an isolated 1 subpattern ---, 0, 0, 1, 0, 0,---. This ef waveform con- 
tains small voltage steps that coincide with the logic level shifts in FE), 
E>, and E3. These steps are not shown in the idealized e} response in 
Fig. 5, but they can occur in a real circuit that has nonzero switch 
offset voltage or stray capacitance feedthrough of switch control logic 
into the signal path. Because the detector bases decisions on the 
difference C;, = B; — Bz-: each cycle, the effects of switch offsets and 
feedthroughs can be eliminated by equating the contributions from 
SWI and SW3 while trimming thresholds I’, and I during 1c manu- 
facture to account for SW2 offset. 
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Fig. 6—Signal and control waveforms for differential detector. 
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Figure 6c shows the peak catcher response e2 after amplification by 
As for the {b,} subpattern ---, 0, 0, 1, 0, 0, ---. The trilevel response 
e3 in Fig. 6b corresponds to this e2 pattern, and is obtained by clamping 
é2 at SW2 when E; is HI. 

Figure 6d gives another comparison of the detector signals for the 
{b,} sequence ---, 0, 0, 1, 0, 0, ---. The test patterns in Fig. 6 were 
obtained for a bubble memory device where no-bubble cycles were 
interleaved with the data cycles. The interleaved no-bubble responses 
are present in the e: waveforms but have been suppressed in ej, és, 
and e3. 


V. CONCLUSIONS 


We have described a magnetic bubble memory detector that will 
make correct decisions from memory sensor outputs even when large 
environmental or chip-to-chip signal transition variations occur in the 
detection interval (D1). The only constraints for successful detection 
are that for the selected functional B, the maximum differences in C; 
= B, — B,-1 for like symbols in contiguous data cycles are less than 
the minimum differences for unlike symbols. When the D1 is properly 
selected to minimize interference, these constraints appear to be ac- 
ceptable for typical outputs from a MBM design when the environmen- 
tal changes in memory operation are slow relative to the field drive 
period. 

At memory start up, the detector must have at least one cycle of 
output from the chip sensors in order to initialize. When input data 
{a,} is encoded to eliminate error propagation, then the recorded 
sequence is {6, }, where b_; = 0, and a single no-bubble sensor response 
is needed to initialize the detector. The comparison of sensor output 
to the uninitialized discharged reference condition in the detector 
circuit in the first read cycle can produce a single error if the first data 
bit in the {6, } data block is no-bubble. This is possible since a memory 
chip may have a no-bubble response B; in excess of the optimum 
differential thresholds [', and I set in the difference detector. This 
detector-initializing no-bubble cycle is not a significant memory system 
constraint. 

The detector circuit considered here incorporates peak sensing in a 
detection interval of the dc restored sensor response to generate the 
B, each data cycle. Other forms of signal processing can be used to 
obtain B,. In particular, if input noise is a problem in future MBM 
designs, the sensor outputs could be dc restored and integrated or 
matched filtered without dc restoration over the detection interval. 
But the resulting “1” and “0” responses for B;, must be sufficiently 
different to guarantee a detectable EYE. Also, if memory device design 
results in excessive intersymbol interference (IsI), signal equalization 
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can be considered to improve detection. At present, peak detection of 
bubble transitions with threshold comparison of C; = B; — B,-1 seems 
to be adequate for differential detection of field-accessed MBM outputs. 
ISI equalization and special care in noise filtering have not been 
required to obtain a consistent differential detection probability of 
error below 10° for the memory types considered. 

In light of the disadvantages frequently encountered with adaptive 
detectors and the advantages of a bubble memory detector that can 
accommodate a wide dynamic range for chip sensor response transi- 
tions, the differential detector described here has considerable merit. 
It is reasonable to integrate this detector as a monolithic bipolar Ic 
such that the bubble memory detection step is accomplished by a 
single self-contained inexpensive DIP. Size and cost as well as the 
detector performance are very important memory system considera- 
tions, since even a modest bubble memory system can require many 
independent detector circuits. 
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Subsampling of a DPCM Speech Channel 
to Provide Two ‘‘Self-Contained’’ 
Half-Rate Channels 


By N. S. JAYANT 


(Manuscript received November 12, 1980) 


We consider the following channel-splitting problem: it is required 
to split a B-bits/s speech-code sequence into two “self-contained” 
B/2-bits/s components, either of which can be used to reproduce 
acceptable speech; also, if both components are available at a re- 
ceiver, it must be possible to reproduce speech with full B-bits/s 
quality. We propose a solution where for interesting values of B, the 
speech quality resulting from half-rate receptions approximately 
equals that from conventional full-rate receptions at B/2 bits/s. In 
the proposed solution, 3.2-kRHz speech is sampled at 12 kHz and 
coded using PcM or differential pcm. The output sequence of code- 
words is split into odd- and even-word sequences. A full-rate receiver 
with access to both of the subchannels simply reconstitutes the output 
sequence prior to decoding, while a half-rate receiver with only the 
odd (or even) subchannel estimates the even (or odd) components by 
nearest-neighbor interpolation. 


1. INTRODUCTION 


The channel-splitting problem described in the abstract is redefined 
in Fig. 1. The receiving end of a speech communication system is 
supposed to operate in either a full-rate or half-rate mode depending 
on whether it has available to it both or only one of the speech sub- 
channels. Respective qualities of speech reproduction are denoted by 
Q(B) and Qu(B/2). The nonavailability of one subchannel is a good 
model for certain types of transmission failure, examples of which are 
signal fading in mobile radio and speech segment losses in packet 
switching. With appropriate forms of diversity reception, the second 
subchannel will be available with probability close to unity when the 
first subchannel is not. The channel-splitting problem has been re- 
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Fig. 1—Definition of the channel-splitting problem. When both of two half-rate 
components in the transmitted sequence are available at a receiver, conventional full- 
rate receptions result, with speech quality Qr(B). When only one of the components is 
present, a half-rate receiver recovers an approximation to the full-rate speech, with 


quality Qu(B/2). 


cently analyzed for communication systems operating at rate-distor- 
tion limits with binary and Gaussian input signals.’” 

The nontrivial nature of the channel-splitting problem can be ap- 
preciated through the simple example of a uniform quantizer. By 
combining two appropriately staggered R-bit quantizers (one of them 
a midrise, the other a midtread), one can realize an (R + 1)-bit system, 
but not a 2R-bit system. For full-rate speech quality corresponding to 
8-bit quantization, component quantizers would each need 7-bit (not 
8/2 = 4 bit) resolution for the combination to yield 8-bit quality. Thus, 
if the subchannels in Fig. 1 were simply uniform quantizers, and if 
speech were sampled at 8 kHz, one would need two 56-kbit/s quantizer 
systems so that a full-rate combination with 64-kbit/s quality can be 
realized. By contrast, in the differential pulse code modulation (DPCM) 
system proposed in this paper, the component receivers that combine 
to give 64-kbit/s quality are indeed half-bit rate, 32-kbit/s systems. 
Moreover, with an illustrative sentence-length speech input, the half- 
rate quality Q:i(B/2) will be shown to exceed the full-rate quality 
Qr(B/2) of a conventional DPcM system operating at B/2 bit/s, for 
interesting values of B. The quality Qr (32 kb/s) is quite acceptable 
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for speech communications, although somewhat short of toll quality. 
Note once again that two uniform quantizers operating at 32-kbit/s 
(8 kHz x 4 bits) each can only give, in combination, 40-kbit/s (8 kHz 
x 5 bits) quality, and not the desired 64-kbit/s quality. Similar argu- 
ments apply to lower bit rates as well. 

The system of Fig. 1 is a special case of a communication scenario — 
that can be generalized to include more than two subchannels, and/or 
subchannels that are non-equal-rate. The system of Fig. 1 can also be 
regarded as a special, symmetrical case of embedded coding,”** with a 
hierarchy consisting of two equally significant subcodes, viz., the half- 
rate sequances. 


Il. SUBSAMPLING AND INTERPOLATION 


The utility of subsampling and interpolation has been demonstrated 
recently in the context of speech packet losses;> speech-encoder out- 
puts are partitioned into odd-sample and even-sample systems which 
are transmitted as separate packets. In the event of a lost odd (or 
even) packet, the lost samples are estimated using nearest-neighbor 
interpolations involving available even (or odd) samples. With the 
usual assumption of 3.2-kHz speech and 8-kHz sampling, the 1:2 
subsampling (at 4 kHz) implies serious aliasing effects, but these errors 
are mitigated by an adaptive interpolation procedure where nearest- 
neighbor-weighting coefficients are varied to follow speech statistics, 
as reflected by appropriate extra information in packet headers.’ The 
system realizes dramatic improvements with packet loss probabilities 
up to about 10 percent; but as the component message loss probability 
approaches 100 percent, as in the channel-splitting problem, residual 
aliasing effects are quite unacceptable, even with adaptive interpola- 
tion. 

The above observation has led us to the notion of 12-kHz sampling 
for the problem of Fig. 1. With 1:2 subsampling, the half-sampling rate 
will now be 6-kHz, which turns out to be just adequate for the 3.2-kHz 
speech inputs in telephony. We also considered 16-kHz sampling, but 
this is less preferable from the point of view of quantization noise. In 
64-kbit/s decoding for example, 4-bit quantization of 16-kHz speech 
produces more quantization noise than the 5-bit quantizer that is 
possible with 12-kHz speech. 

A second advantage of 12-kHz sampling is that it permits nonadap- 
tive interpolations; adaptive interpolations yielded near-zero gains 
with 12-kHz speech. If 8-kHz is subsampled, it cannot be adequately 
reconstituted by nonadaptive interpolation even if the interpolation is 
invoked with a probability much less than 100 percent.° 

Waveform reconstructions in the half-bit-rate (half-sample rate) 
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systems of this paper are described by 
u(n) = A-u(n — 1) + Az-u(n + 1), 
A; = Az = 0.5. (1) 


The samples u(r) will be quantized DPcM prediction error samples 
in general; and more simply, they will be quantized speech samples in 
the special case of nonpredictive, or nondifferential Pcm. 


lll, THE DPCM CODEC 


Figure 2 shows block diagrams of full-rate and half-rate DPCM 
systems with fixed first-order predictors. In each case, decoding is 
defined by 


y(n) |ppom = hi-y(n -l1)+ q(n), (2) 


where qg(n) is the quantized prediction error signal, and h; is a first- 
order predictor. In the special case of nondifferential Pcm (h; = 0), 
q(n) is simply the quantized speech output: 


y(n) |rom = q(n). (3) 


Subscripts H and F in Fig. 2 distinguish half- and full-rate versions of 
y(n) and q(n). 
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Fig. 2—prcm block diagrams. (a) Conventional full-rate codec. (b) Decoder portion 
of half-rate receiver. 
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In contrast with well-known adaptive DPCM systems where the 
quantizer step size is adapted for every sample, the present paper 
assumes a system where the step size is adapted once for each block 
(of several ms duration), and held fixed for the duration of the block.® 
A periodically updated, rather than instantaneously adaptive, quan- 
tizer is used in anticipation of interpolation procedures, which are 
known from recent experience to be unreliable when the step size 
exhibits sample-to-sample fluctuation.’ 

The periodically adaptive quantizers are defined by a block-specific 
step size A that is proportional to the root-mean-square value of a 
generalized first-difference in the form 


A |aqr = Ki(R)-[x(n) — hi-x(n — 1) rms, (4a) 
A lags = Ki(R)-[y(n) — Ai-y(n — 1) ] ms, (4b) 

with maximum and minimum constraints 
Amax = 128-Amin = Ko(R)-|x(7) |max; (4c) 


where [ Ki(R), Ko(R)] are bit-rate-specific constraints with suggested 
values of [0.25, 0.03], [0.33, 0.06], [0.58, 0.12] and [1.0, 0.18] for 5, 4, 3 
and 2-bit quantizers respectively. The subscripts AQF (4a) and AQB 
(4b) refer to forward-adaptive and backward-adaptive procédures; 
respective rms values in (4a) and (4b) are evaluated over the duration 
of a speech block to be coded in AQF, and over the duration of the 
most recent decoded speech block in AQB. The AQB procedure is less 
effective because of speech nonstationarity as well as the effect of 
quantization noise that is present in the y(n) sequence used in (4b). 
However, step-size information in an AQB system need not be sepa- 
rately transmitted to a receiver; it is inherently available in the decoded 
y(n) sequence. AQF procedures, by constrast, require the explicit 
transmission of step-size information (typically, about 5 bits worth, per 
block of 16 ms). In our experiments, quality losses in AQB were more 
noticeable in full-bit-rate speech than in half-bit-rate speech; and in 
each case the losses were of a second order of importance. With this in 
mind we have elected to cite only AQF results in section IV; these 
results can be regarded as upper bounds as far as quantizer perform- 
ance is concerned. 

In the context of 1:2 subsampling, the AQB procedure of (4b) cannot 
be implemented as such unless h; = 0 (PcM). However, step sizes 
obtained by setting fh; to zero in (4) have béen found to have fairly 
small effects on ppcM performance. Differences between PcM-optimal 
and DPcM-optimal step size are less significant than differences among 
step sizes of different speech blocks. Moreover, the suboptimality of 
PCM-matched step size becomes less significant as h; — 0. The next 
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section shows that half-bit-rate Dpcm favor values of h; that are indeed 
much smaller than those appropriate for conventional full-rate pPcM. 


IV. RESULTS AND CONCLUSIONS 


Figures 3 through 6 illustrate the performance of the interpolation 
procedure for the example of a 5-bit encoder. The waveform segments 
refer to two 20-ms blocks from a 3.2-kHz bandlimited, female utterance 
“The chairman cast three votes” sampled at 12 kHz. 

Figure 3 shows full-rate and interpolated gq waveforms (fA; = 0.5) for 
the two segments; it demonstrates that the nonadaptive interpolator 
(1) is reasonably adequate even for the fast-varying unvoiced example. 
This is confirmed in Fig. 4 which shows corresponding full-rate and 
half-rate y waveforms. Notice that the half-bit-rate output provides a 
much better reproduction of the voiced segment, but is nevertheless 
reasonably effective in unvoiced speech reproduction. Perceptually, 





(b) 


Fig. 3—Full-rate and half-rate (interpolated) sequences qr and qi (R = 60 kb/s, Ai 
= 0.5) for (a) voiced speech segment and (b) unvoiced speech segment. The speech 
segments are 16 ms long. 


506 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1981 





(b) 


Fig. 4—Original and half-rate-decoded speech output y(n) (R = 60 kb/s, A: = 0.5) for 
(a) voiced speech segment and (b) unvoiced speech segment. The speech.segments are 
16 ms long. 


the waveform degradations in the half-bit coder are indeed fairly subtle 
for R = 5. 

The significance of h; = 0.5 is demonstrated in Fig. 5, which shows 
DPCM performance as a function of predictor coefficient value. The 
objective quality measure used is the segmental signal-to-noise ratio 
SEGSNR defined as the average of 10 log s/n (dB)-values measured over 
the totality of 20-ms blocks in the input. Notice that maximization of 
full-rate and half-rate quality call for h; = 0.9 and hi; = 0.5, respectively; 
and notice also that these are not very sharp maxima, suggesting 
flexibility for practical implementations. The special, simple case of hi 
= 0 (pcm) leads to a noticeable quality degradation only for the full- 
rate system. 

Figure 6 depicts the performance of full-rate and half-rate DPCM 
receivers as a function of encoder bit rate. The performance curves (1) 
and (ii) are for h, values that maximize half-rate speech quality (/1 
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Fig. 5—Segmental sNR versus predictor coefficient h; (R = 60 kb/s). Values that 
maximize full-rate and half-rate quality are respectively 0.9 and 0.5. 


= 0.5 for 5,4 bits/sample, h, = 0.3 for 3,2 bits/sample). The full-rate 
characteristic shows the expected 6-dB-per-bit behavior, while the 
half-bit-rate characteristic falls more gradually with decreasing R. 
Both characteristics tend to the expected 0-dB limit for no transmis- 
sions (& — 0). The square dots in Fig. 6 represent the performance of 
a full-rate receiver in a system designed to maximize full-rate speech 
quality (h; = 0.8 to 0.9). 

An important observation from Fig. 6 is that for encoder bit rates in 
the important range of 2 to 5 bits/sample, 


Qu(R/2) = Qr(R/2); R < 5 bits/sample. (5) 
This suggests that the half-bit-rate qualities in the subsample-inter- 
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refer to half-rate and full-rate codecs designed to maximize half-rate quality. The square 
dots represent a full-rate codec designed to maximize full-rate quality. 
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polate system are extremely good results, considering the crucial 
constraint that the half-bit-rate systems combine trivially to yield the 
full-rate performance Q(R). The approximate equality in (5) is borne 
out very well in perceptual assessments of Qy and Qr. In contrast to 
(5), analytical results in Refs. 1 and 2 are quite pessimistic. This 
difference in conclusions is related to the fact that these analytical 
results apply at the rate-distortion limit, while the bit rates in this 
paper are nowhere close to the rate-distortion limit for speech. In fact, 
our bit rates are high enough that there is sufficient redundancy left in 
the coder output to permit subsampling and high-quality interpola- 
tions. 

The relative performance of the half-bit-rate receiver diminishes 
with increasing bit rate. Clearly, as R — ©, the quantization noise 
contributions to @ vanish, Qr— ©, and Qy tends to a finite asymptotic 
value that shows the effect of nearest-neighbor interpolation noise. 
Results elsewhere’ can be used to show that this asymptotic value, for 
a first-order Markov signal example, is approximately given by the 
expected value of 10 log [(1 + R2x(1))/(1 — R2(1)] dB, where Rx (1) 
is a block-specific adjacent sample correlation in the speech input 
x(n). 

Finally it would be appropriate to calibrate the Qr and Qu values in 
Fig. 6 with well-known definitions of toll-quality and communications 
quality (near perfect intelligibility with noticeable but not obstrusive 
degradation).® Although simulations and conclusions have centered on 
a single earlier-cited test input, it appears that full-bit-rate DPCM 
realizes toll-quality for R = 5 and 4 bits/sample; and communications 
quality for R = 3 and 2 bits/sample. The half rate DPcM receptions in 
the proposed system approach toll quality with R = 5 bits/sample and 
maintain good communications quality at R = 4 and 3 bits/sample. 
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Blocking States in Connecting Networks 
Made of Square Switches Arranged in Stages 


By V. E. BENES 
(Manuscript received October 16, 1980) 


Since the probability of blocking is a principal measure of the 
performance of a network, this paper examines the blocking states of 
a network. For two-sided connecting networks made of square 
switches arranged in stages, a parallel pair (PP) is a pair of paths 
through the network that meet at no switch. The blocking states of 
the network are closely related to those that have a busy pp. In 
particular, a routing algorithm can avoid all the blocking states if 
and only if it avoids all the states with busy PPs. 


l. INTRODUCTION 


Connecting networks made of rectangular switches have long been 
an important feature of telephone central offices, and their combina- 
torial properties still command substantial practical as well as theo- 
retical interest. The probability of blocking, calculated or measured in 
a suitable setting, is a principal measure of the performance of a 
network, and it is of particular interest to identify and study those 
combinatorial features which give rise to blocked calls. 

To be sure, blocking is known to be an unfortunate mismatch of 
available resources. But precisely, what structural features lead to it? 
What presages blocking? What are some tell-tale signs and necessary 
concomitants of blocking? When and how can blocking be avoided? 
To try to answer these questions, we focus on the blocking states of 
the network, those in which some call is blocked. 

Although in a century of telephony much work has gone into 
calculations and simulations for blocking probability, there are (except 
for a few examples) virtually no results on the basic causes and nature 
of blocking. For networks made of square switches arranged in stages, 
this paper will partly remedy this lack by introducing some simple 
concepts suggested by examples; theorems based on these bear on the 
position and structure of the blocking states within the set (semilattice) 
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of states, and also on the possibility of routing rules that avoid blocking 
entirely. 


il. PRELIMINARIES 


In a network made of square switches each rectangular switch has 
the same number of inlets as outlets, so there is neither expansion nor 
concentration. The network is said to be arranged in stages if the 
switches are partitioned into sets called stages, two of which carry the 
inlets and outlets, while the others are internal stages, such that every 
path from an inlet to an outlet passes through each stage only once, 
always in the same order. The network is full access if each switch in 
a given stage has a link to every switch in the adjacent stages. 

It has been pointed out repeatedly’”® that the set S of permitted 
states of a connecting network forms a partially ordered system, in 
fact a semilattice, in which the order relation x <= y between states x, 
y means that x can be reached from y by zero or more hangups: In this 
case we say that x is below y. The usual Hasse diagram for the partial 
ordering is then a'convenient way of visualizing the set S of states, 
marred only by the prodigious number of states that arise. This “curse 
of dimensionality” can be alleviated by considering only equivalence 
classes of states under a suitable group of symmetries,” and we do so 
when convenient without further discussion. 

A state x € S is maximal if there are no states above x. It can be 
seen that in a network made of square switches all maximal states 
carry the same number of calls. A state x is a blocking state (or simply, 
is blocking) if an inlet u and an outlet v exist, both idle in x, such that 
there is no available allowed path from u to v, whose use would connect 
u to vu and give rise to a state y covering (immediately above) x in the 
partial ordering; the call (u, v) is then said to be blocked in x. We use 
the convention that, in speaking of a call (u, v), u is always an inlet, 
and v always an outlet. Figure 1 illustrates all the preceding concepts. 

A routing algorithm or rule is a method for choosing routes for 
unblocked calls during the operation of the network to carry traffic. It 
specifies, for each state x and each call (u, v) not blocked in x, what 
route should be used to put up (u, v) if these terminals desired 
connection when the system was in state x. Needless to say, it is of 
interest to find and describe rules which make for good network 
performance as evaluated by the probability that an arriving call is 
blocked, i.e., the fraction of requests for service which cannot be 
completed because of blocking. Since such denials of service occur only 
when the network is in a blocking state, it is natural that a good 
routing strategy will seek to avoid the blocking states. Examples are 
known in which the blocking states can be avoided entirely; such 
networks are called nonblocking in the wide sense.”” 
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INLET STAGE OUTLET STAGE 





Fig. 2—A parallel pair (PP). 


Ill. SUMMARY OF RESULTS 


We assume as given a two-sided connecting network made of square 
switches arranged in stages, and we introduce the idea of a parallel 
pair (PP). This is a pair of paths through the network that meet at no 
switch (Figs. 2, 3). The basic message of this paper is that visiting 
states with busy pPs is bad for performance, because the property of 
being a blocking state is closely tied up with that of having a busy PP. 
Now having a blocked call and having a busy PP are usually not 
equivalent properties of states; the characterization of blocking is not 
that simple. There are nonblocking states which have a busy PP, and 
there are blocking states without a busy pp. Nevertheless, our intuitive 
conclusion is strongly supported by small examples,** and here by the 
following results: 

(t) Every maximal state that is above a blocking state has a busy PP. 
(tz) Every blocking state is below a maximal state with a busy PP. 
(it) A routing algorithm makes the network nonblocking in the wide 
sense if and only if it makes inaccessible all states with busy PPs. 

It can be seen that in a state without any busy PPs, every path in use 
must meet every other at a switch. To maintain this property in an 
operating system, one must keep from adding paths that pass only 
through empty switches. Thus, the advice to avoid busy Pps almost 
coincides with (or at least, is very consistent with) the “packing” 
principle that to avoid blocking one should route calls in the most 





Fig. 3—Paths that meet at the second stage. 
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heavily loaded part of the network that will accept them: Clearly, in 
putting up a new call you best avoid busy Pps by having the route you 
choose meet as many calls in progress as possible. Our result says that 
completely avoiding blocking and never having any busy PPs are both 
possible, or both impossible, to achieve by a choice of routing rule. 


IV. SOME EASY RESULTS 


To fix ideas and record a few facts, we first list various results with 
straightforward proofs. The reader interested in warming up to the 
concepts in use here might go through these in some detail using Figs. 
1 to 5. 

(1) Remark: No state immediately below a maximal state is block- 
ing. 

(2) Remark: A blocking state two levels down from a maximal state 
(i.e., obtainable from a maximal state by two hangups) consists of an 
idle pp with all other terminals and links busy. 

(8) Lemma: If x is a blocking state, then some PP is idle in x. 

Proof: Let (u, v) be blocked in x. Since the switches are square, 
there is an idle path from wu to the other side, and one from vu to w’s 
side of the network. These paths must avoid each other if (u, v) is 
blocked. Thus, they constitute a pp idle in x. 

(4) Remark: If x has a busy pp, then states y and z exist such that 
y is maximal, z is blocking, y => x, and y = z. 

Proof: Starting in state x, put calls in until a maximal y is reached, 
then hang up the PP to form z, and reverse its assignment, 1.e., if the 
busy PP provided the two calls (uw, v;) and (ue, U2), then the calls (uw, 
v2) and (us, v,) are both blocked in z. Compare Remark (2). 

(5) Remark: If v is full access and x is blocking, then all inner 
switches are nonempty in x. 

Proof: An empty switch is accessible from any switch in the adjacent 
stages, hence from any idle inlet, and also from any idle outlet. Thus, 
no call is blocked. See Fig. 4. 

(6) Remark: x nonblocking + x has no busy PPs. 

Proof: Figure 4 provides a counterexample. 
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Fig. 4—State with two busy Pps but no blocking. 
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(7) Remark: x blocking + x has a busy PP. 
Proof: Figure 5 gives a counterexample. 


V. BLOCKING STATES ARE BELOW STATES WITH BUSY PPS 


(8) Theorem: A maximal state with no busy Pps has no blocking 
states under it. Every maximal state that is above a blocking state 
has a busy PP. 

Proof: We embroider the argument in (3). Suppose that x is a 
blocking state, and y = x is maximal. Let (u, v) be an inlet-outlet pair 
blocked in x; u and v are on disjoint paths of a pp idle in x. Moreover, 
since y = x, it is possible to start with x and put in all the additional 
calls of y, by the same routes as they use in y, except those involving 
the terminal pair (u, v). (This is not to say that a given routing 
algorithm would prescribe those routes in operation!) Since (u, v) is 
blocked in x, wu and v are not talking to each other in y. Indeed, in y, u 
is connected to some switch other than u’s, and v to some switch other 
than w’s. If the paths carrying these calls in y met at some switch, (u, 
v) would not have been blocked in x. Hence, in the state formed from 
x by putting in all the calls of y except those involving u and uv, these 
paths again form an idle pp with wu on one path of the PP, and v on the 
other. Thus, y necessarily has a busy PP. 

(9) Remark: We have already seen in (4) that a maximal state with 
a busy PP is above some blocking state. Thus, for maximal states y, 


y has a busy PP iff y lies above a blocking state. 


(10) Corollary: Any blocking state lies below a maximal state 
which has a busy PP. 

Proof: Every state x lies below some maximal y. If y had no busy 
pps, no blocking state could be below it, by (8). 

(11) Corollary: If x is below some maximal state which has no busy 
pps, then x is nonblockhing. 

Proof: (8) or (9); the property is sufficient, but not necessary. 


VI. THE CLOSURE TOPOLOGIES INDUCED BY = 


A partial ordering <= induces a topology in a familiar way. If X C S 
is a subset of states, its lower closure, denoted Xic, is the set of states 
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Fig. 5—Blocking state with no busy PPs: (u, v) is blocked. 
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that are below some member of X in the sense of S, including, of 
course, the elements of X as well. Lower closure satisfies the Kura- 
towski axioms for a topological closure operation. There is a dual 
topology of upper closure, wherein the dual = replaces <=, denoted by 
aes 

By introducing notation for some of the important subsets of states 
we have encountered so far, we can use lower closure to express various 
preceding results in a concise way. Let us define 


Max = set of states that are maximal in <, 
B = set of blocking states, 
S — B = set of nonblocking states, 
bPP = set of states that have a busy parallel pair. 


(12) Remark: (bPP),, = (Max N bPP)t. 

Proof. If x € (OPP), then there exists y with a busy pp such that x 
< y. Pick a maximal z = y. Then z € OPP, since bPP is closed above. 
Thus, we have z = x and z € Max N OPP, and so also x € (Max N 
bPP),-. The converse is trivial. 

In these notations, Theorem (8) may be cast as saying that 


(13) BN (Max N (S — bPP),. = § = empty set, 
while Corollary (10) claims that 
BC (Max N bPP)., 


which by (12) simplifies to B C (bPP))-. Similarly Corollary (11) can 
be rendered as 


[Max N (S — bPP))}. CS — B. 


Vil. WIDE-SENSE NONBLOCKING NETWORKS 


A telephone network is called nonblocking in the wide sense (NBWs) 
if a rule exists for assigning routes to incoming calls that precludes all 
blocking, that is, under the rule no blocking state is accessible. Run 
according to the rule, the system can satisfy all demands for connection 
as they arise without ever rejecting any, or rerouting existing calls. 
The point is to try to reduce network complexity (from the extreme of 
strict nonblocking) by using a network that has blocking states, but to 
avoid them by clever routing. 

Two necessary and sufficient conditions are known for a general 
connecting network to be NBws: One is phrased in terms of denseness 
in a simple metric on the states,” the other in terms of a semilattice 
homomorphism of S onto the set of assignments realized by the states 
in S.° For networks made of square switches arranged in stages, we 
give a new necessary and sufficient condition, based, perhaps not 
surprisingly now, on avoiding states with busy PPs. 

If a is a routing algorithm, we use Ac(a) to denote the set of states 
that are accessible under a, i.e., that can be reached by the operating 
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network if it is run using a to make routing choices. A good routing 
algorithm derives its worth from making bad states unlikely or inac- 
cessible, so it is important to know what states can be reached. Further, 
it is clear than an NBWS network is one for which an a can be found 
such that Ac(a) N B is empty, i.e., no blocking state can be reached by 
operating under a. We note these properties of Ac(a): 

(14) Ac(a) is closed below: [Ac(a) ]ic = Ac(a), 

(15) x € Ac(a) = Ac(a) N Max N {x} #8. 

Property (15) says that if x is accessible, then some maximal state 
above x is accessible. For networks made of square switches, however, 
a specialized sharp form of (15) holds, in which both hypothesis and 
conclusion are stronger, and which relates blocking to having a busy 
pp. To be specific, we shall use a precise version of the following 
heuristic idea: An algorithm a must complete an unblocked call by 
some route; so starting in a state x that has a blocked call (u, v), it 
should be possible to keep picking successive unblocked new calls not 
involving u or v, and routing them by a, until as many terminals are 
busy as possible, and the system has arrived at a state z accessible 
from x under a by new calls, with (u, v) idle and blocked in z; we claim 
in particular that the calls can be picked so that z consists of an idle 
PP which has u and vu on switches that terminate its paths, and all 
other terminals are busy. Compare Remark (2). 

(16) Lemma: If (u, v) is a call blocked in state x, then an inlet s 
and an outlet t exist, with s ¥ u and t ¥ v such that EITHER 

(i) -s, tu, and v are the only terminals idle in x, and there is a PP 
idle in x one of whose paths joins s to v, the other u to t, OR 

(ii) an inlet q and an outlet r exist, neither of which coincides with 
any of s, t, u, or v, such that q and r are both idle in x and (q, r) is not 
blocked. 

Proof: We recall our convention that in a call (u, v), uw is always an 
inlet and v an outlet. Since (u, v) is blocked, there is a PP idle in x one 
of whose paths connects u’s switch to some switch other than v’s, 
which, since the switches are square, has an idle outlet ¢ on it. Similarly, 
the other path of the pp connects v’s switch to some switch other than 
u’s on which there is an idle inlet s. Suppose that alternative (z) fails, 
so that some inlet qg, distinct from each of s and u, is idle in x. Now 
since switches are square, g can reach an idle outlet r; we shall show 
that r can be chosen to be distinct from both ¢ and v. There are three 
cases. 

Case (i): uw and g are on the same switch. Since uw can reach t, q 
must be able to reach some r other than ¢. This is because if we put up 
(u, t) in x, g could still reach some idle outlet, so a fortiori it can do so 
in x. Since (u, v) is blocked, so is (q, v). But (qg, 7) is not blocked, so r 
~ v. Thus, r is distinct from both ¢ and v, though it may be on the 
same switch as t. So alternative (it) holds in Case (1). 
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Case (ii): s and qg are on the same switch. Then g can reach v’s 
switch and possibly ?¢’s. If it can reach v’s switch only, then an idle 
outlet 7 must exist on u’s switch, r ¥ v. This is because if an (s, v) call 
were put up, g could still reach some idle outlet, which must be on vu’s 
switch, since adding a call does not make additional outlet switches 
reachable by qg. Thus, r ¥ t¢ and r ¥ uv, (q, r) is not blocked, and 
alternative (ii) holds. Suppose then that g can reach ?’s switch as well 
as u’s. There are two subcases: (a) ¢’s switch has another idle outlet r 
distinct from ¢, or (b) ¢ is the only idle terminal on ?¢’s switch. In case 
(a), r#v, rt, and (q, r) is not blocked, so we are done. In case (5), 
since both wu and s can reach f¢, the final arc in the paths from u to ¢ 
and s to ¢ must be the same. Thus, these paths meet at an inner switch, 
and indeed on the way through the stages from inlets to outlets they 
meet for the last time at an inner switch o. Thus, o has at least two 
idle terminals on each side, and there is an idle path from o to an outlet 
switch other than ?’s, which is also reachable by u, because the paths 
meet. This outlet switch is therefore distinct from vu’s because (u, v) is 
blocked; it has an idle outlet 7 on it distinct from each of ¢ and v, and 
(q, r) is not blocked. Thus, again alternative (zi) holds. 

Case (iii): g is on neither w’s switch nor s’s. If g can reach any switch 
other than ?¢’s or vu’s, we are done; for then there is an idle outlet r, 
distinct from each of ¢ and v, and (q, r) is not blocked. The same is 
true if g can reach only ?’s or u’s switch provided the one it can reach 
has more than one idle outlet. Suppose then that qg can reach only ?’s 
switch or u’s, and that whichever ones of these it can reach have 
exactly one idle outlet, viz, ¢ or v as the case may be. 

If g can reach ¢ but not v, then (q, v) is blocked and moreover there 
is an inner switch at which q’s idle path to ¢ meets one from u to ¢ for 
the last time (indeed they coincide henceforth), and there is a path, 
subsequently disjoint from the remaining common path to ¢, from that 
inner switch to some outlet switch other than ?’s or v’s, because (gq, v) 
is blocked. On that outlet switch is the r we seek for alternative (ii). If 
q can reach v but not t, the same conclusion follows, with s playing the 
role of u, mutatis mutandis. 

If g can reach both ¢ and uv, we argue similarly, as follows: Since ?’s 
switch has exactly one idle link, an idle path from u to ¢t and one from 
q to ¢ must have that last idle link in common, there being no other 
way to reach ?’s switch. Hence, these paths must meet at an inner 
switch at least once. As before there is then a last, in the ordering of 
stages from the inlet side to the outlet side, switch at which these 
paths meet, and they coincide thereafter. Since the paths are meeting, 
the switch has at least two idle links incoming from the inlet side, and 
at least two idle links outgoing on the outlet side. One of the latter is 
a part of the now common part of the paths on their remaining legs to 
t’s switch. The other idle link on the outlet side is part of an idle path 
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to the outlets; it does not meet the common part again, and so cannot 
end at ?’s switch. If it reached vu’s, (u, v) would not be blocked. The 
outlet switch it reaches then has an idle outlet r on it, distinct from 
each of ¢ and v, and (gq, r) is not blocked; hence alternative (ii). 

(17) Lemma: If x is a blocking state accessible under a routing 
algorithm a, then some maximal state y = x is accessible under a, 
and y has a busy pp. In symbols, 


x€ BN Ac(a) > Ay y=x, y € Max N Ac(a) N OPP. 


Proof: Let (u, v) be a call blocked in an accessible state x. If possible, 
use alternative (ii) of Lemma (16) to find a call (q, r), unblocked in x, 
neither of whose terminals is on an idle PP with u and v. Putting up 
(qg, r) according to a results in a state in which (u, v) is blocked. Add 
calls in this way repeatedly until no such (gq, r) can be found, i.e., 
alternative (ii) fails. Then alternative (i) holds: The state z we have 
reached consists of an idle pp with one of u, v on each of its distinct 
paths. It is obvious that there is one maximal y above z, reachable 
under a by having the right two calls request connection, and that y 
has a busy PP. 

(18) Theorem: Let a be a routing algorithm. Some blocking state 
is accessible under a iff some state with a busy PP is accessible under 
a; in symbols, 


Acla)N B#®H iff Acla)N bPP AQ. 


Proof: Suppose x € Ac(a) N bPP. By (15) there exists y => x such 
that y is maximal and y € Ac(a). Since y = x and bPP is upper closed, 
we see that y has a busy PP. By (9) there is a z € B such that z <= y and 
z € Ac(a), since Ac(a) is closed below. Hence, Ac(a) N B ¥ ¥. The 
converse follows from (17). 

We illustrate two aspects of Theorem (18) by reference to Figs. 4 
and 5. According to the converse part of (18), if an algorithm a made 
the state in Fig. 4 reachable, then there must be a continuation (of the 
events that led to this state) leading to a blocking state. In Fig. 4 this 
happens if there is, first, a call between the two highest outer switches, 
and then, a call between the two second-highest outer switches. There 
are two ways of satisfying this extended assignment, and both lead to 
blocking. The direct part of (18) is exemplified in Fig. 5. If an algorithm 
a makes the state in Fig. 5 reachable, then according to (18) a state 
with a busy PP is also reachable. It can be verified that in Fig. 5 putting 
up any new call leads at once to a state with a busy PP. 

(19) Corollary: A routing algorithm a makes the network nonblock- 
ing in the wide sense if and only if no state with a busy PP is accessible 
under a. 

An example in which it actually happens that a natural routing rule 
makes all blocking states unreachable was shown in Fig. 1. The relation 
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B shown in Fig. 1 by arrows is to be interpreted as “better than.” Any 
algorithm for routing which uses the “better” alternatives avoids all 
blocking and all busy pps. Note that there are actually two such 
algorithms, as suggested by the “don’t care” label. 
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Optical Fibers for Scanning Digitizers 
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Automated graphic activities such as computer cartography and 
computer-aided design require document-scanning systems that are 
stable, fast, and capable of handling large drawings and maps. In 
such a system the scanning assembly should accept documents readily 
and be compact and inexpensive enough to encourage wide usage. 
An attractive method for such automatic scanning and digitizing of 
engineering drawings combines optical fibers with a solid-state im- 
aging array and a stepped friction feed. In this paper we describe an 
assembly of these elements employing a light-guide illuminator and 
receiver that are in direct contact with the drawing. Light signals, 
modulated by the dark areas on the drawing, are transmitted from 
the linear end of the lightguide to a grid configuration at the opposite 
end. The signals are optically coupled to and detected by a solid- 
state array sensor that scans the grid of optical fibers. The dimen- 
sional control necessary for fabricating the fiber grid to match the 
configuration of the solid-state array is critically important for the 
operation of this type of scanner. We therefore emphasize in this 
paper the dimensional requirements for the glass fiber assembly, the 
signals detection, and the accuracy of scanning. 


|. INTRODUCTION 


Document scanning systems for automated activities such as elec- 
tronic mail, computer cartography, and computer graphics require an 
image-guiding system to capture visual details for conversion to elec- 
tronic signals. The image arrangement is combined closely with, and 
influences the design of, the document-positioning mechanism. The 
document can be illuminated by diffused light in the general scanning 
area, the bounced light optically projected through a tiny aperture 
onto an appropriate sensor (flying aperature scanners), or by spot 
lighting where a tiny spot of light is projected onto the copy and the 
bounced spot is picked up directly by the sensor (flying spot scanners). 
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The surface condition alters the light that is reflected to the photoe- 
lectric detector during the scanning. Documents are scanned electron- 
ically, mechanically, or in combination, with the portion of the docu- 
ment being scanned mounted flat or in a circular arc. 

The physical dimensions of scanning assemblies are controlled by 
the size of the document and the light paths required to scan the 
surface. In most assemblies, lenses and mirrors are combined with 
rotating or translating document holders.’ In more recent assemblies, 
stationary copy is scanned by a flying spot from a cathode-ray tube,’ 
a charge-coupled photodiode array,” or an acousto-optically modulated 
laser.** Optical fibers are also used to pipe light from the document to 
a rotary scan head or solid state image sensor.’ In this paper, we 
describe the performance characteristics of a flying-spot type of optical 
fiber scanner that has inherent advantages for digitizing large docu- 
ments. This paper gives special attention to glass fiber characteristics 
required to construct this new potentially low-cost scanning digitizer. 


ll. OPTICAL FIBER SCANNER 


A digitizer must scan the document to resolve the line and character 
information into small picture cells or pixels. The optical response at 
each pixel is then used to control detection circuitry to generate 
positionally defined digital output. In general, the drawing can be 
scanned in one of three ways: individual pixels can be addressed in a 
series of parallel sweeps across the drawing surface; groups of pixels 
can be addressed by a sensor array that is stepped over the drawing; 
pixels can be addressed sequentially in parallel sweeps to locate lines 
and marks on the drawing that are subsequently followed or tracked. 
In all methods, some form of parallel scanning must first be made to 
locate the line information to a prescribed resolution. As will be shown 
for scanning large documents, optical fibers, suitably arranged for a 
raster output (method 1), offer high resolution, geometric fidelity, 
stability, and mechanical simplicity—all in a very compact assembly. 

Figure 1 shows the schematic representation of an optical fiber 
scanning digitizer. Light is transmitted from the source via a multi- 
plicity of lightguides to uniformly illuminate a narrow strip across the 
drawing. The fibers at the light source are randomized in a compact 
circular group and are illuminated with a small lamp. The receiver 
consists of a similar number of optical fibers as the illuminator. These 
are oriented in a line at one end and as a matrix or grid at the other 
end. The linear end of the illuminator and receiver are geometrically 
similar and located adjacent to each other in the fixed scanner head. 
Drawings are positioned directly against the fixed head by a spring- 
loaded pressure plate. The drawing is moved past the read head in 
small increments by rubber-coated friction rollers driven by a preci- 
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Fig. 1—Schematic representation of an optical fiber scanning digitizer. 
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sion-stepping motor. The stepping motor, operating at a few millisec- 
onds per incremental movement, is driven by a sensor controller. The 
paper is moved when the controller signals that it has completed a 
linear scan of the document. The diffused and reflected light signals, 
attenuated by dark lines, characters, etc., on the drawing, are trans- 
mitted through the lightguide receiver and detected by an array of 
photosensitive elements that scan the lightguides. A one-to-one cor- 
respondence exists between the specific transverse location on the 
drawing, a receiver optical fiber, and an element in the solid state 
array. 

A light mask and lenses are interposed between the receiver grid 
and the image sensor array. The light mask, prepared by a photolith- 
ographic process, is a metalized glass plate with an equal number of 1- 
mil-diameter holes as receiver fibers. The light-mask holes are located 
on the nominal centers of the fibers. The mask corrects for slight 
misalignments in the fiber grid and presents a dimensionally precise 
pattern to the array sensor. The lenses focus the light from each fiber 
(through the light-mask holes) onto the individual sensor elements 
(pixels). Gray scale or thresholded video output from each pixel is used 
to represent the digitized signals. These signals can be processed 
further® to provide information that is suitable for reproduction, com- 
puter editing,’ or graphics teleconferencing.® 

Table I gives a list of the important characteristics of an optical fiber 
scanning digitizer. Figure 2 shows the assembly used to scan digitize 
by optical fibers, which provided the experimental data described in 
this paper. 


itl, FIBER ASSEMBLY DESIGN 


In scanning digitizers utilizing optical fibers, fabricating the fiber 
grid so that its tolerances match the location of the active pixels in 
solid state array sensors is critically important. While in communica- 
tion applications fiber diameter is an important consideration regard- 
ing transmission parameters, in scanning digitizers by optical fibers it 
is crucial that the fiber grid, hole mask, and solid state array are 
correctly aligned. The degree of precision that can be achieved in 


Table |—Characteristics of an optical fiber scanning digitizer 


Type Optical fiber/ccp array/roll feed 

Copy Width Up to 40 in. 

Copy Length Unlimited 

Copy Feed Variable, 1 ms/0.004 in. maximum 

Scan Resolution Variable, dependent on fiber diameter and spacing 
Maximum Pixel Rate 10 mH (present maximum) 

Output Gray scale or threshold video 

Size Desk top (10 in. high, 20 in. deep, and 50 in. wide) 
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positioning the fibers in a grid depends on the dimensional control of 
the fiber elements and the mode of fabrication of the assembly. 

Problems arise in the construction of scanner grid assemblies that 
are similar to problems with fiber optic array connectors. In two- 
dimensional array connectors, the alignment mechanism is either a 
hole mask through which the fibers are threaded or thin chips are used 
that are grooved on both sides and stacked to hold the fibers.? These 
methods work well and are economical for the relatively small number 
of fibers in cable. However, for aligning fibers for scanners that use 
many thousands of fibers, construction processes should be automated 
to avoid handling individual fibers. A method of forming plastic optical 
waveguide by selective photopolymerization” offers one means of mass 
fabricating large linear arrays.’ Unfortunately, the plastic sheets 
would cause fabrication problems. For example, dimensional variabil- 
ity of the sheet thickness would cause an imprecise alignment of the 
lightguides in the grid. A method that can be used successfully involves 
wrapping continuous fiber on a drum with a transverse groove for the 
subsequent injection of adhesive to fuse the assembly to form precise 
linear arrays. Further, the linear array can be split into bands of equal 
amounts of fibers for subsequent stacking into grids of the required 
configuration. 

Figure 3 shows dimensional requirements for individual fibers and 
for the cross section of bands of fibers that are assembled by the drum 
method. In this example, there are 100 fibers in each band (250 fibers/ 
in.). The 0.004 in. diameter (101.6 pm) is typical for optical fibers 
currently fabricated for communication systems, but the required 
tolerance of +0.00001 in. (+0.25 ym) standard deviation is well beyond 
normal practice. To achieve this precision requires using advanced 
fiber drawing control” and instrumentation systems that use noncon- 
tacting on-line diameter measuring.’*"* As shown in Fig. 3, fibers in 
the band are uncoated and are positioned in tangential contact with 
each other. Epoxy fills the interstitial spaces on the bottom of the 
band. The adhesive joins the fibers together and forms a flat surface 
that is used to support the stack of other bands in making the grid. 

Figure 4 shows an enlarged portion of a grid assembly. In this 
assembly, the vertical dimension is also controlled by the fiber diam- 
eter. A second application of adhesive, this time to fill the remaining 
interstitial spaces between fibers, is used to secure the bands to form 
a compact and precise grid. As noted previously, a light mask placed 
over the grid can compensate for minor misalignment induced by slight 
variations in fiber diameter and processing. The light mask consists of 
an array of 10,000 0.001-in. holes, located on 0.004-in. centers. The 
mask is fabricated by the same technology used for solid state sensor 
arrays. The final assembly therefore provides a grid of tiny emitting 
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Fig. 3—Dimensional requirements for optical fiber linearly aligned lightguides at 250 
bits/inch resolution. 


holes with dimensional control similar to that obtained in solid state 
array fabrication. 


IV. SCANNER PERFORMANCE OF FIBER OPTICS 


Figure 5 shows the results of scanning a test pattern to evaluate the 
characteristics of a fiber assembly made by the drum and band method. 
Lines on the test pattern, varying in width from 0.0040 in. to 0.0145 in. 
produce detectable signals via the 0.004-in. diameter fibers. Fibers 
used in this test are of stepped-index type with a 0.7 core to cladding 
diameter ratio. The fibers are spaced nominally on 0.004-in. centers. 
As shown in the figure, this dimensional arrangement of fibers can 
resolve into discrete signal/line widths that are well below the 0.015 
in. minimum that is the ANSI Y14.2 reconimendation for engineering 
drawings. 
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Fig. 4—Portion of optical fiber grid receiver. 


In Figure 5 the digitized signals represent a partial scan across the 
drawing in the X direction. The absolute accuracy of scanned data in 
the X direction (perpendicular to the document motion), relative to 
the graphics on the drawing, is controlled by the positional alignment 
of the optical fibers as shown in Fig. 3. In the experimental unit each 
fiber is positioned within +0.0005 in. to adjacent fibers, and each band 
of 100 fibers varies in width by at most +0.0005 in. However, the 
placement of bands adjacent to each other to form the linear portion 
of the lightguide receiver introduces additional variations that cause 
apparent stretching or shrinking in the width of the engineering 
drawing. This effect can be compensated by an appropriate step in 
subsequent data processing to define the physical location of each 
lightguide rather than to assume the nomimal spacing. 

Accuracy in the Y-direction scan is controlled by the rotational 
characteristics of the stepping motor and dimensional variation in the 
friction feed rollers. Commercially available gear-driven stepping mo- 
tors operate with an accuracy of +4 arc min. This results in a maximum 
single-step error of 0.00026 in. per step, but because of gearing, the 
error is noncumulative. Dimensional variations measured in experi- 
mental friction-feed rollers can contribute less than 0.2 percent varia- 
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tion to the individual step distance. This can produce a cumulative 
error and results in an apparent stretching or shrinking of the drawing 
in the Y direction. To obtain a more accurate digital rendition of the 
drawing, this condition can also be corrected by data processing. 
Figure 6 shows the output from an experimental scanning digitizer 
for a portion of an engineering drawing. This display of thresholded 
video, as recorded in sections on a cathode-ray tube, shows the digiti- 


ACTUAL WIOTHS OF LINES AT 
CROSS SECTION 





SCANNED 
L “CROSS SECTION 


i THRESHOLDED VIDEO ‘Ss 


aeons OOOCOOHHEHRHREERERS FESEESHE RoR eee EE Fee OEE RED eet EtereEe eaves 





Fig. 5—Signals from a 250 fiber/inch scanner. 
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Fig. 6—Digitized output from an experimental optical fiber scanner as recorded in 
sections on a cathode-ray tube. 


zation in the X and Y direction. Irregularities in the boundaries of the 
lines of the drawing appear as steps because of the digitization process. 
Files of this type can be processed further by computer for a variety of 
automated graphics activities. One use of this type of output is to 
create graphic files for computer-aided design systems.” 


V. CONCLUSION 


This paper has presented an overview of the considerations and 
requirements in the design of optical fibers for the scanning and 
digitizing of large documents. It shows that when combined with a 
solid state array and a roller feed, this type of assembly has distinct 
advantages. High resolution (250 bits per inch) can be obtained with 
fiber that has a diameter that is typical for communication systems. 
Geometric fidelity is inherent in the digitized output of this type of 
scanner because of the fixed-dimensional arrangement of fibers and 
the stepping motor-driven rollers. These aspects also contribute to the 
stability and mechanical simplicity of the assembly. Because the fibers 
are flexible, they can be folded into a compact arrangement to produce 
a desktop-size assembly for scanning the widest of engineering draw- 
ings. This type of assembly also permits convenient handling of the 
drawing with entry and exit from the front of the scanner. 

Of critical importance to optical fiber scanning digitizers is the 
fabrication of the fiber grid to tolerances that match those of solid 
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state array sensors. For very large grids construction processes must 
be automated to avoid handling of individual fibers. We have found 
that wrapping continuous fiber on a drum with a transverse groove for 
the subsequent injection of adhesive forms precise and stable linear 
arrays. Further, the linear array can be split into bands of equal 
amounts of fibers for stacking and fusing into a grid. To achieve the 
required precision for the fibers, instruments that use noncontacting 
on-line diameter measuring, as well as advanced drawing and feedback 
controls, are necessary. 

Using optical fibers for scanning digitizers can open new possibilities 
for computer-aided design and engineering. The availability of a tech- 
nology that permits compact, high-resolution, and stable scanners 
should encourage wider use than now exists for automatic digitizing. 
Coupled with advanced data compaction software, convenient-to-use 
low-cost scanners can serve aS a new input means for graphics pro- 
cessing systems and can permit the conversion to data of the largest of 
existing and newly prepared engineering drawings. 
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A magnetometer array has been installed aboard the SCARAB 
submersible to enable it to detect and track buried undersea telephone 
cable. Four three-axis magnetometers sense the magnetic field of a 
25-Hz current applied to the cable. The magnetometer signals are 
filtered, amplified, digitized, and multiplexed onto SCARAB’s umbillt- 
cal cable for transmission to the control ship. A shipborne minicom- 
puter processes these signals in real time to determine the cable 
location, which is continuously displayed to SCARAB’s operator on a 
graphics terminal. This paper describes the development and capa- 
bilities of the system. Among the topics discussed are a history of the 
cable-locating problem, an analysis of the factors governing the 
achievable 25-Hz signal level, a description of the magnetic noise 
spectrum and its sources, and a discussion of the signal-processing 
techniques. We also examine the dependence of system performance 
on signal and noise levels. 


1. INTRODUCTION 


SCARAB (Submersible Craft Assisting Repair and Burial) is a new 
class of unmanned cable-controlled submersible craft whose mission is 
to facilitate the maintenance of undersea telephone cable systems. The 
first SCARAB, now under development, is owned and will be operated 
by a consortium including AT&T Long Lines and British, Canadian, 
and French telecommunications companies.* A second SCARAB belongs 
to Transpacific Communications, Inc., a subsidiary of AT&T Long 
Lines. The SCARAB vehicle, shown in Fig. 1, is typical of other recently 
designed unmanned submersibles. The basic structural element is a 


* Transpacific Communications, Inc. (a subsidiary of AT&T Long Lines); Cable and 
Wireless, Ltd.; French Postes Télégraphe et Téléphonie; British Post Office; Teleglobe 
Canada. 
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Fig. 1—SCARAB’s major subsystems. 


frame of welded aluminum rod. Flotation tanks are mounted on the 
frame to provide buoyancy, and seven thrusters are used for propul- 
sion. Two of the thrusters are oriented vertically to allow the vehicle 
to dive and ascend. Two manipulators are provided for cutting and 
gripping cable. A water pump is used in conjunction with dredger and 
jetter attachments for unburying and burying cable. Additional equip- 
ment includes TV systems, acoustic devices, and a sensor system that 
includes the cable detection apparatus. 

Operationally, SCARAB is designed to perform two specific functions: 
to assist during the repair of faulted cable, and to perform surveys of 
operational cables. During a repair operation, SCARAB magnetically 
locates the cable track and, subsequently, the point of the fault. It then 
unburies and cuts the cable as required, and attaches grippers so the 
cable ends can be lifted to the cable ship. Once both ends are aboard 
the ship, the cable is spliced, tested, and then lowered to the sea bed. 
The scARAB then cuts away the lowering lines and reburies the exposed 
length of cable. 

Basically, the survey operation is to follow the cable and measure 
the burial depth with the use of the magnetometer detection system. 
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If a point of shallow burial or cable exposure were detected, the vehicle 
would rebury using the jetter. 

The purpose of the present paper is to describe in some detail 
SCARAB’s cable-locating subsystem, designed and constructed by Bell 
Laboratories. More precisely, SCARAB has two distinct cable-locating 
systems, for reasons discussed below. The ac system described herein 
is considered the primary system; the dc system (not described) is 
intended to serve as a backup under certain operational circumstances. 


ll. HISTORY OF CABLE-LOCATING PROBLEM 


Telephone cable systems use single coaxial cables whose outer 
diameters are between 2.5 and 5 cm, with repeaters inserted every 5 to 
40 nautical miles, depending on the system. (The newer systems have 
the larger cable and closer repeater spacing.) Since the first transat- 
lantic telephone cable was laid in 1956, the major threat to undersea 
cables has been the trawling and dredging operations of fishermen on 
the Atlantic continental shelves. In fact, about 95 percent of undersea 
telephone cable failures have been attributed to fishing activities. The 
principal response to this threat, aside from more heavily armoring 
the shelf sections of cables, was the development in 1966 of the first of 
a series of plows (Sea Plows I to IV) for burial of ocean cable.’ Early 
versions of the sea plow buried the cable to a depth of 70 cm, leaving 
the repeaters exposed; the latest plow” buries the repeaters as well. 

Burial reduces the likelihood of cable faults but by no means 
forestalls them. Even a buried cable system is vulnerable near unburied 
repeaters, or where cable burial was impeded by adverse bottom 
conditions, or where bottom erosion occurs after plowing, or where 
ships’ anchors may drag. 

This vulnerability, though minimized, nevertheless poses a problem: 
if a buried cable is damaged, how can it be found and its fault 
localized? An approximation to the fault location can be obtained from 
shore by combining knowledge of the cable’s route with an estimate of 
the offshore distance to the fault, the latter obtained by any of several 
electronic fault-locating techniques applied at the cable’s shore ter- 
minal. This procedure can identify the faulty repeater section, thereby 
isolating the fault to within a few miles, but the problem remains for 
the cable-repair ship to pinpoint the fault location and recover the 
cable. This problem is considerably exacerbated by burial of the cable, 
a fact recognized at the time the first sea plow became operational. To 
provide repair ships with the capability of finding buried cables, Bell 
Laboratories developed three varieties of electromagnetic probe: an 
electrode pair to be streamed behind the repair ship with a horizontal 
spacing of about 50 m; a 76-cm long magnetic induction coil to be 
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towed at low speeds in a vertical orientation near the sea floor; and a 
smaller induction coil for hand-held use by a diver. 

To use the probes, a low-frequency signal (20 to 40 Hz) is applied 
between the cable’s center conductor and ground at the shore cable 
terminal nearer the fault. This signal propagates along the cable, and 
its electric field can be detected at sea by the electrode system 
streamed behind the repair ship, allowing the ship to track the cable 
at relatively high speeds by zigzagging across it. Unless signal atten- 
uation is excessive, the most seaward location at which the signal is 
detected provides a crude fault-location estimate, which can be refined 
by towing the induction coil at slow speeds so that it hangs nearly 
vertically, indicating by a null when it is above the cable. It was 
intended at this stage of the repair operation that a diver, using the 
hand-held probe to pinpoint the fault location, attach lifting lines and, 
if necessary, cut the cable. However, because of the depth limitations 
and inherent danger of diving operations, the usual procedure is to 
pick up the cable by grappling. Although grappling is quite straight- 
forward for unburied cable, grappling has obvious disadvantages when 
applied to buried cable. The grapnel may not penetrate the bottom 
consistently to the cable’s burial depth, thereby running the risk of 
missing the cable completely; or the cable may not have adequate 
tensile strength to withstand being pulled obliquely (often nearly 
horizontally) out of the bottom through more than a meter of mud or 
sand. These problems and similar ones motivated SCARAB’s develop- 
ment. 


lil. CABLE LOCATING BY SCARAB: OPERATIONAL REQUIREMENTS AND 
SYSTEM CONCEPTS 


3.1 SCARAB’s cable-locating requirements 


During a typical repair operation, SCARAB is called upon to find a 
buried cable and track it for possibly a mile or more until the cable 
fault is located. Thereupon it must position itself on the bottom, close 
to the cable, so its dredger and manipulators can be used to unbury 
the cable, cut it if necessary, and attach a lifting line. Cable reburial 
follows cable repair. SCARAB may also be required to estimate the 
cable’s burial depth, both during repair operations and during survey 
of operational cable, to determine whether the cable should be buried 
more deeply. 

In view of the complexity of these tasks and the difficulty of remotely 
maneuvering a swimming vehicle within a few meters of the bottom, 
the cable-locating system aboard SCARAB must be considerably more 
sophisticated than the ship-towed system described above. A principal 
limitation of the ship-towed system, from the viewpoint of its possible 
application aboard SCARAB (e.g., by mounting a single induction coil 
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on the vehicle), is its reliance on a single sensor, which measures only 
one component of the electric or magnetic field vector, and at only one 
point in space. Clearly, one instantaneous measurement of this kind is 
insufficient to localize the cable, which therefore can only be found by 
maneuvering the probe to a variety of positions near the cable (e.g., 
zigzagging), noting the field component as a function of time (and 
therefore of space), and using this time series to estimate the cable’s 
location. 

This procedure serves well where the information required concern- 
ing the cable’s location is minimal and qualitative, such as during a 
grappling operation, where one usually needs only to know which of 
two opposite courses will cross the cable. However, for quantitative 
cable localization (e.g., to determine cable burial depth), this procedure 
would of course require real-time computer processing; but further- 
more the processing algorithm would require knowledge of the vehicle’s 
track over the bottom so that past history could be used to determine 
present cable location. Hence the computer would need accurate data 
concerning vehicle speed and ocean current; of these, the former is 
difficult to measure accurately at slow speeds and the latter is virtually 
impossible to measure from a swimming vehicle. (The problem is 
analogous to ship navigation using running fixes, a procedure whereby 
past and present lines of position are conjoined to determine the ship’s 
position. The past line of position is updated using a dead-reckoned 
ship’s track, which is subject to errors in the estimated course and 
speed of the ship and the set and drift of the current. The procedure 
is viable only because the ship’s speed is usually large compared with 
the errors in its assumed speed over the bottom, and because the 
requirements for positional accuracy are not stringent. Neither of these 
helpful conditions pertains to SCARAB’s cable-localization problem.) 

So that estimates of track over the bottom do not have to be relied 
upon, a precision cable-locating system must have enough sensors so 
that their outputs, taken at a single instant of time, are sufficient to 
determine cable location without using any past data. For this purpose, 
a minimum of five sensors is required, since four parameters are 
needed to specify a given line in space (the cable location) and one 
parameter is needed to specify the cable current. In fact, for reasons 
discussed below, it proves convenient to use 12 sensors, packaged in 
four three-axis modules. 


3.2 Choice of sensor 


An attractive possibility for a cable-locating sensor is a bottom- 
penetrating sonar. In fact, a parametric sonar has succeeded in detect- 
ing telephone cable buried under 30 cm of sand in a laboratory tank 
experiment.” Although the feasibility of a sonar cable detector has 
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been thereby strongly suggested, the time and funding required to 
develop an operational sonar system for SCARAB were judged to exceed 
the available resources. Accordingly, attention was focussed on passive. 
electromagnetic sensors. 

The next choice to be made was between electric and magnetic 
sensors. Because at least five sensors would have to be mounted on a 
small vehicle crowded with other equipment, compactness was an 
overriding concern. From this viewpoint, electric sensors (i.e., electrode 
pairs) are a disadvantage because their size, determined by their 
baseline length, is normally considerably larger than a comparably 
sensitive magnetic detector. Furthermore, the technology of commer- 
cially available magnetic sensors is well advanced, having benefited 
from long-standing geophysical and military interest in compact and 
reliable magnetic detectors (an interest with no counterpart in electric- 
field detection). For these reasons, and because the geometry of the 
cable’s magnetic field offers some processing advantages over its 
electric-field geometry, magnetic cable detection was chosen. 

Magnetic detectors fall into three categories, according to the mag- 
netic quantity measured: total-field magnetometers (measuring the 
field magnitude), component magnetometers (measuring individual 
vector field components), and induction coils (measuring the time 
derivative of individual field components). 

Of these, induction coils can be readily dismissed from further 
consideration. Their sole advantage over magnetometers is their lower 
self-noise level, which is irrelevant in this application because other 
noise sources dominate self noise. Their principal disadvantages are 
the size and the weight of their cores (important considerations for a 
small equipment-laden vehicle that must be positively buoyant). 

The choice between total-field and component magnetometers is 
not as straightforward. To clarify the functional difference between 
them, consider the case of a total-field magnetometer in the presence 
of the earth’s field B. and a cable field B., where | B.| « | B.| (as is 
always the case in practice). If b is a unit vector parallel to B., then 
the magnetometer output is | B. + B.| = | B.| + b-B., which is simply 
the component of the total field parallel to the earth’s field (i.e., the 
perpendicular components of the cable field do not contribute). A 
total-field magnetometer is therefore completely equivalent to a com- 
ponent magnetometer oriented along the earth’s field. 

Suppose now that B, is time-varying, so that its contribution to the 
magnetometer output can be separated from that of the earth’s static 
field by filtering, leaving a signal of b-B. from the total-field device 
and a-B. from the component device (where 4 is a unit vector parallel 
to the magnetometer axis). The sole distinction between these outputs 
is that b is fixed while 4 is changeable (i.e., a component magnetometer 
can be oriented to measure any desired component). 
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This flexibility of the component magnetometer is both an advantage 
and a disadvantage. The advantage is that three component magne- 
tometers with mutually orthogonal axes can be combined in a single 
package to create a vector sensor (i.e., a sensor which detects all three 
components of the vector field). Because the electromagnetic field is 
a vector field, there are crucial signal-processing advantages to be 
gained by employing a vector sensor, as will be seen. The disadvantage 
of the component magnetometer is that it is subject to a potentially 
severe source of noise if mounted on a moving vehicle. As the vehicle 
maneuvers, pitches, and rolls, it changes the orientation of the mag- 
netometer’s axis relative to the earth’s field, so that in the earth 
coordinate system a is time-varying. The magnetometer output is 
therefore [B.(¢) + B.]-a(¢), and the contribution from the earth’s field 
can no longer be removed by filtering unless a(t) has negligible 
frequency components in the frequency band of the signal B,(¢). 
Because the earth’s field normally exceeds the cable field by nearly 
four orders of magnitude, even slight vehicle motions can produce 
intolerable levels of motional noise if these motions are in the signal 
frequency band. 

The choice between total-field and component magnetometers 
therefore hinges on whether a usable signal frequency exists, at which 
motional noise is acceptably small. As will be shown below, SCARAB 
has sufficient inertia and stability, so when it is submerged, motional 
noise is negligible at frequencies above a few Hertz. This fact was not 
known when the choice of sensor had to be made, because SCARAB had 
not yet been built; accordingly, an experiment was conducted using 
the Navy’s CURV III, a cable-controlled underwater recovery vehicle 
quite similar in design to SCARAB. A component magnetometer 
mounted aboard CURV III revealed not only that motional noise was 
acceptably low at all frequencies of possible interest, but also that 
magnetic noise from all sources (notably the thrusters) was low enough 
to encourage the development of the cable-locating system described 
below. 


3.3 Choice between ac and dc signals 


A further preliminary question concerned whether the system should 
operate at ac or dc. An ac system would detect the ac magnetic field 
of a low-frequency current on the cable. A dc system would detect the 
superposition of dc fields from two sources: dc current on the cable, 
and permanent magnetization of the cable’s steel components (central 
strength member or external armor wires). 

The overriding difference between ac and dc cable detection derives 
from the existence of a large dc component of the earth’s field. To 
detect the small dc cable field in the presence of the earth’s field, one 
must use the fact that the earth’s field is spatially homogeneous, 
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whereas the cable field is a localized anomaly. Therefore, to detect the 
cable one must detect differences between field values at different 
locations in space; if these values differ, the cable is nearby. This 
procedure can be implemented with one sensor by moving it in space 
and watching for changes in signal. Or it can be implemented with 
multiple sensors by looking at the differences among their signals. The 
first of these techniques must employ a total-field sensor, since a vector 
sensor would be subject to changes in orientation as it was moved (the 
dc equivalent of motional noise). The second of these techniques could 
employ either total-field or component sensors, though usually com- 
ponent sensors are used, with two being packaged in the same housing 
to maintain their parallelism. Such a configuration is termed a gradi- 
ometer, since its output is proportional to the field gradient if the 
sensor separation is small. . 

One contrast between ac and dc cable detection, therefore, is in the 
type of sensor that can be used: ac detection can use a vector magne- 
tometer, while dc detection requires a total-field magnetometer or 
gradiometer. Put another way, the ac sensor measures a vector quan- 
tity while the dc sensor measures either a scalar or a tensor quantity. 
Because the magnetic field is a vector, one might expect the vector 
sensor to lend itself more readily to data processing for cable localiza- 
tion, thereby favoring ac over dc. In fact, this expectation is borne out: 
the ac cable-locating algorithm described below is extremely simple, 
whereas a comparable dc algorithm has not yet been devised. 

Nonetheless, a dc cable-locating system possesses certain advantages 
over ac. Since undersea telephone cables normally carry a dc current 
(with sea return) to power the repeaters, a dc detection system does 
not require a special signal to be placed on the cable, and the system 
works as well for an operational cable as for a broken cable. Further- 
more, unlike an ac signal, the dc current does not diminish exponen- 
tially with distance from shore, so both ends of a broken cable can be 
powered with comparable current levels. In fact, when the cable is 
armored it does not need to be powered at all for dc detection, which 
can take advantage of the relatively large fields of the steel armor 
wires. In contrast, an ac signal cannot propagate across the ocean, and 
therefore, at best, it can reach only one side of a break in a transoceanic 
cable. Under certain circumstances, therefore, the unique capabilities 
of a dc system may be necessary. 

In sum, dec sensors provide magnetic data of a type that is not 
suitable for cable localization, with the result that an ac system is 
normally preferable. However, operational circumstances exist for 
which an ac system would be ineffective and a backup system is 
desirable. For these reasons, a dc system is under development for 
SCARAB to supplement the ac system. 
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The remainder of this paper describes the implementation and 
operational performance of the ac system. 


3.4 Cable-localization algorithm and array design 


In the preceding sections we proposed that the optimum cable- 
locating system should employ a vector sensor (i.e., a sensor detecting 
all three field components), because in some sense a vector sensor is 
better suited to characterizing a vector field. A corollary was that the 
optimum system was an ac system, to avoid the earth’s dc field. A 
further proposition was that the optimum system should have enough 
sensors to determine cable location instantaneously, without having to 
rely on past history. 

To particularize these general considerations, this section describes 
the method implemented aboard scARAB for instantaneous cable lo- 
calization by a number of ac vector sensors. The method relies on the 
fact that the magnetic field lines caused by the cable’s current are 
circles (assuming the cable lies in a straight line). This geometry 
permits cable localization by a simple triangulation procedure most 
easily visualized in two dimensions (Fig. 2). Vector sensors are used to 
determine the direction of the field at two locations. At both locations 
a line perpendicular to the field is constructed. Because these two lines 
lie along radii of concentric circles, their point of intersection is the 
common center: the cable location. The generalization of this proce- 
dure to three dimensions is straightforward: at two locations, construct 
planes perpendicular to the field direction, and determine the line at 
which these planes intersect. This line of intersection is the cable’s 
location. 

In principle, this procedure allows cable localization with only two 
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vector magnetometers. Moreover, even this basic system overdeter- 
mines the problem: a line in space can be described with four param- 
eters, whereas two vector sensors provide six scalar measurements. 
These six measurements can be manipulated to provide information 
on field direction (four angles) and size (two magnitudes). Cable 
localization requires only the four angles; the two magnitudes can 
provide two cable-current estimates, which would of course coincide 
in the absence of noise and other sources of field distortion. 

However, this simple system of two vector sensors must be embel- 
lished, for it fails when the two sensors and the cable all lie in the 
same plane (or, in the presence of noise, nearly in the same plane). 
This degenerate case cannot be ignored, for it occurs quite commonly 
for any two-sensor array. For example, to optimize cable tracking it 
would be logical to place the two sensors symmetrically to port and 
starboard, so that the sensors straddle the cable when SCARAB does. 
For such a placement, however, SCARAB would be unable to do cable 
localization while heading perpendicular to the cable, such as during 
initial cable acquisition. 

The solution is to provide more than two sensors, and to implement 
a procedure for deciding which sensors shall be used and which ignored. 
To ensure that no cable location shall be coplanar with all pairs of 
sensors, a nonplanar array of at least four sensors is required. 

The optimum position of these four sensors is governed by the desire 
that the planes whose intersection determines cable location should 
intersect nearly at right angles, since then the cable-location errors 
caused by noise-induced fluctuations of these planes will be minimized. 
For right-angle intersection, the baseline length between two sensors 
must be twice the distance from the baseline to the cable. However, 
SCARAB’S maximum baseline length is constrained to be less than two 
meters (since booms extending beyond the vehicle’s perimeter are 
undesirable), while SCARAB’s height above the sea floor will generally 
be at least one meter. Consequently, right-angle triangulation is im- 
possible to achieve, and the effects of noise are minimized simply by 
spacing the sensors as far apart as possible, maximizing the distance 
between each sensor and the plane defined by the other three. In fact, 
even this modest goal is impractical, because most potential sensor 
sites in the center and after parts of the vehicle are contaminated by 
thruster noise and frame-current distortion of the magnetic field (see 
below). Consequently, the best available array configuration comprises 
four coplanar sensors in a nearly vertical plane in the forepart of the 
vehicle (Fig. 1). This array is not well suited to localizing a cable which 
lies athwartships and underfoot (i.e., a cable in the plane of the 
sensors), but for all other cable locations it serves well. 

Such an array of four sensors provides six possible sensor pairs and 
therefore six cable-location estimates, if the outputs from all pairs are 
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used. Since some pairs will be better situated than others to triangulate 
the cable, a choice must be made among them. This choice is made by 
calculating a weighted average of the six cable-location estimates, 
which de-emphasizes sensor pairs whose measured magnetic fields are 
nearly parallel. The weighting function used is sin 0,* where @ is the 
angle between the magnetic fields and therefore between the planes 
whose line of intersection is sought. This weighting function satisfies 
the requirement that it increase monotonically from 0 at @ = 0° (the 
most undesirable case) to 1 at @ = 90° (the most desirable case), but 
beyond this requirement its choice is rather arbitrary. 

This section has described how SCARAB’s ac cable-location system 
works, in principle. Outputs from an array of four vector magnetome- 
ters are combined in a shipboard computer, using a rather simple 
cable-locating algorithm and weighting scheme, and the resulting 
estimate of cable location is displayed to SCARAB’s operator so that he 
may maneuver the vehicle accordingly. How well this sytem works in 
practice depends on the relative levels of signals and noise. This is 
addressed in Section IV. 


IV. SIGNAL AND NOISE 
4.1 Signal 


The magnetic field H of a line current J in an infinite homogeneous 
nonconducting medium satisfies 


|H}=—, (1 
2p 

where H is in amperes/meter, J is in amperes, and p is the distance 
from the line current in meters. In a conducting medium such as sea 
water, this expression is valid only if p < 6, where 6 is the skin depth. 
In sea water, however, 5 = 250/f'”” meters, where fis the frequency in 
hertz, so that 6 > 40 m for all frequencies less than 40 Hz. Equation 
(1) is therefore well satisfied within about 20 m of the cable. This 
conclusion is unaffected by the presence of the sea bottom, whose 
lower conductivity causes no appreciable field distortions except at 
distances greater than a sea-water skin depth. An important implica- 
tion is that the circular geometry of the field lines, on which SCARAB’s 
cable-location method relies, is not appreciably distorted by the pres- 
ence of the sea bottom except at ranges in excess of 20 m. 

To characterize the signal, therefore, it remains only to specify J, 
the net current on the cable. When a signal current is fed to the center 
conductor at the shore terminal, return currents flow partially in the 
sea and partially in the outer conductor, metallic shielding tapes, and 
armor wires, if any. The net current J, which determines the field 
according to eq. (1), is the algebraic sum of all currents on the cable’s 
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various conductors. This net current attenuates as the signal propa- 
gates seaward, because of losses in the cable, the sea water, and the 
repeaters (which introduce loss rather than gain because the low- 
frequency signal follows the dc power path through the repeaters 
rather than the path through the carrier-frequency amplifiers). For 
example, Fig. 3 shows a theoretical prediction of propagation at several 
frequencies along a 100-mile section of armored cable with 5-mile 
repeater spacing and with a cable break at the offshore end. The 
appendix describes the theory underlying Fig. 3. The net current is 
lower at higher frequencies, not only because the losses are greater, 
but also because the outer conductor carries a greater share of the 
return current (so that the net current, being the algebraic sum of 
center- and outer-conductor currents, is less). Clearly, from the view- 
point of signal propagation, we prefer the lower frequencies. 

The choice of operating frequency, however, is ideally based on 
maximizing the signal-to-noise ratio. This ideal requires knowledge of 
the vehicle’s noise spectrum, however, and SCARAB’s development 
schedule did not permit delaying the choice of frequency until the 
vehicle was completed and its noise spectrum measured. Therefore, 
we provisionally chose 25 Hz for convenience, since 25 Hz is a frequency 
commonly used with ship-towed probes (though its choice in that 
context was influenced by considerations entirely irrelevant to SCARAB, 
such as the low-frequency falloff of induction-coil sensitivity and the 
expected low-frequency noise caused by turbulence-induced motion of 
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Fig. 3—Propagation of low-frequency signals on armored sc cable with a cable break 


100 nmi offshore. The center-conductor current is 200 mA at the shore point, for all 
frequencies. The step discontinuities occur at the repeaters, which are unpowered. 
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Fig. 4—Propagation at 25 Hz on armored sc cable with cable breaks at various 
locations. The unlabeled curve extending to 220 nmi is for an unbroken cable. The 
center-conductor current at the shore point is 200 mA (but note that the vertical scale 
extends only to 80 mA). The step discontinuities occur at the repeaters, which are 
unpowered. 


a ship-towed coil). Although the choice of 25 Hz was therefore arbitrary 
to some degree, subsequent determination of SCARAB’s noise spectrum 
(Section 4.2) showed that 25 Hz lies in a rather quiet region of the 
spectrum, and so we have retained it as SCARAB’s operating frequency. 

Figure 4 illustrates propagation at 25 Hz on the same type of cable 
system as Fig. 3, but with cable breaks at various distances offshore. 
Notice that a cable break enhances the signal by total reflection of the 
propagating mode, thereby roughly doubling the cable current at the 
fault relative to the current on an intact cable. A broken cable is 
therefore twice as easy to detect as an intact one (if both carry the 
same current at their shore points). 

Aside from the question of signal attenuation, the maximum current 
that can be propagated to the desired detection point is restricted by 
the maximum center-conductor current tolerated by the repeaters. 
This value is 100 or 200 mA (rms), depending on the cable system, 
except when the system is operational, in which case no more than 20 
mA (rms) can be tolerated without affecting service. Broken cables 
therefore can carry at least ten times the signal of operational cables, 
with a factor of at least five due to the higher allowed current and a 
factor of two due to the signal reflected from the cable break. 


4.2 Noise 


Knowledge of the noise field and its sources is a prerequisite to 
predictions of operational performance, and also may suggest methods 
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of performance improvement. SCARAB’s magnetic noise originates from 
a variety of sources, described below in rough order of decreasing 
importance. 


4.2.1 Thruster noise 


SCARAB’s propulsion is by thrusters comprising induction motors 
(2300 volt, four pole, three phase, 60 Hz) driving shrouded screws. Six 
5-hp thrusters are used, four for fore-and-aft propulsion and two for 
vertical propulsion (Fig. 1). A seventh thruster, horizonally oriented 
athwartships, is powered by the hydraulic system, whose pump is 
driven by a 35-hp motor. The noise spectra of CURV’s and SCARAB’s 
thrusters were measured in tank tests and, as would be expected for 
any rotating machinery, consisted of spectral lines with no measurable 
broadband component. 

The basic frequencies of an induction motor’s operation are the 
shaft-rotation frequency and the synchronous frequency (i.e., the 
apparent rotation rate of the stator’s magnetic field). For a thruster 
with no external load or internal friction, these frequencies are the 
same, equalling the power frequency (60 Hz) divided by the number of 
pole pairs (two for SCARAB, three for cuRV). Under load the thruster 
slows down, reducing the shaft-rotation frequency by an amount 
which, when multiplied by the number of pole pairs, is called the slip 
frequency. For example, the synchronous frequency of the SCARAB’s 
thrusters is 30 Hz, while the shaft and slip frequencies at one-third 
speed are 10 Hz and 40 Hz, respectively. 

Thruster noise is principally at four frequencies: the slip and shaft 
frequencies, and their sum and difference frequencies. Noise at higher 
harmonics, although measurable at high speeds, is small. The geometry 
of the noise field at the slip and shaft frequencies is consistent with an 
oscillating dipole source aligned parallel to the thruster axis. In con- 
trast, the geometry of the noise field at the sum and difference 
frequencies is consistent with an oscillating and rotating dipole source 
aligned perpendicular to the thruster axis. By far the strongest of these 
lines is the slip-frequency line. As the speed of SCARAB’s thrusters 
varies between zero and 90 percent speed (the maximum attainable 
under load), the frequency of this line varies between 60 Hz and 6 Hz, 
and its root-mean-square amplitude (expressed as an effective source 
dipole moment) varies between 0 and 19.6 y-m* (where 1 y = 1/400 z 
amp/m). The amplitude is roughly a linear function of the voltage, 
and therefore varies quite nonlinearly with motor speed (at half speed, 
the noise amplitude is about one quarter of its full-speed value). 

The precise origin of these noise components is obscure, related in 
some way to secondary leakage fields and therefore to subtle details of 
motor design and manufacture. An illustration of the subtlety is that 
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a 5-hp SCARAB thruster has about six times the noise of a 10-hp CURV 
thruster when both are operated at maximum speed. 

Thruster noise poses a potentially awkward problem, for it is com- 
posed of several spectral lines whose frequencies can attain any value 
between 0 and more than 60 Hz, depending on the operating speed. If 
the signal frequency were chosen higher than any attainable by thrus- 
ter noise, a heavy penalty would be paid in the attenuation of the 
signal in its propagation from shore (Fig. 3). The alternative is to base 
the choice of frequency on other considerations and rely on being able 
to avoid thruster speeds producing excessive noise at the chosen 
frequency. The strongest constraint is on the use of the forward 
vertical thruster, which is closest to the magnetometer array. The 
operator must therefore be particularly careful in the use of that 
thruster while the cable-localization system is being used. In fact, 
during normal search and tracking operations the operator customarily 
relinquishes control of the vertical thrusters to SCARAB’s automatic 
altitude control system, which servo-controls these thrusters to main- 
tain the desired height above the bottom. Fortunately, this automatic 
altitude control system uses the vertical thrusters in a manner quite 
favorable to the cable-localization system. 

To illustrate the effects of thruster noise, Fig. 5 shows a noise 
spectrum recorded using a three-axis magnetometer mounted on 
SCARAB at a point midway between the two upper magnetometers 
shown in Fig. 1. This point is well within a meter of the forward 
vertical thruster, is about 1.5 m from the 35-hp hydraulic-pump motor, 
and is more than 2 m from the thrusters at the stern. The spectrum 
was recorded during dockside tests in San Diego harbor, while SCARAB 
was swimming submerged at nearly full speed with the vertical thrus- 
ters controlled by the automatic altitude system. We feel that this 
spectrum is reasonably characteristic of the typical noise environment 
of all SCARAB’s magnetometers. 

In Fig. 5, the intense spectral line near 40 Hz comes from the forward 
vertical thruster, whose proximity to the magnetometer causes it to 
dominate the other noise sources. At the time the spectrum was taken, 
this thruster was operating at one-third speed, generating noise at its 
slip frequency near 40 Hz and at its sum frequency (i.e., shaft frequency 
plus slip frequency) near 50 Hz. These values of the frequency reflect 
typical operating speeds of the vertical thrusters, whose average speed 
in level flight is, of course, a direct function of the vehicle’s buoyancy. 
The automatic altitude system maintains level flight using rather small 
speed excursions, in contrast with manual altitude control in which 
the natural tendency of the operator is to overreact by using large 
changes in vertical thrust. 

In comparing the three components of the 40-Hz line, note that the 
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smallness of the athwartships component (relative to the other two 
components) is consistent with the tank-test result that the slip- 
frequency noise source is equivalent to a dipole parallel to the thruster 
axis. 

Other spectral lines in Fig. 5 come from the horizontal thrusters and 
the hydraulic-pump motor. The horizontal thrusters, operating at 
nearly full speed, contribute noise at their slip frequency of about 10 
Hz. Although they are generating more noise than the vertical thrus- 
ters (because of their higher speed and their greater numbers), their 
impact on the noise spectrum is substantially smaller because of their 
greater distance from the magnetometer. 

The 35-hp hydraulic-pump motor operates at a constant speed just 
below its synchronous speed of 30 Hz, and its most important contri- 
bution to the noise spectrum is a strong line at its sum frequency (slip 
plus shaft frequency) just above 30 Hz. Its slip frequency is very low 
(near 1 Hz), and so its slip-frequency noise is of no practical impor- 
tance. 

From the spectrum in Fig. 5 it is clear that 25 Hz is a fortuitous 
choice of operating frequency, since it avoids the frequencies at which 
thruster noise is most severe. The frequencies of the horizontal and 
vertical thrusters are variable, of course, but it is quite practical to 
avoid the critical speed (58 percent) at which their slip frequency is 25 
Hz. The vertical-thruster speed necessary to maintain level flight is 
substantially less than 50 percent at the vehicle’s buoyancy, and the 
automatic altitude system maintains vehicle altitude without ap- 
proaching the critical speed of 58 percent. On the other hand, we can 
expect the speed of the horizontal thrusters to pass through 58 percent 
in normal operation, but their noise level at this speed is substantially 
lower (by 8 dB) than their level in Fig. 5. As long as SCARAB does not 
swim continuously at this speed, the effects of the horizontal thrusters 
should be confined to small, short transients seen by the magnetome- 
ters each time the critical speed is passed. 

In addition to the spectral lines discussed above, the spectrum of 
Fig. 5 contains noise at the power frequency and several other fixed 
frequencies. The 60-Hz line contains contributions from all the vehi- 
cle’s equipment, but the predominant sources are the thrusters and 
the pump motor, since they draw the most current. Subharmonics of 
the power frequency at multiples of 12 Hz are also present, caused by 
unidentified equipment on board the vehicle other than the thrusters 
and pump motor. The subharmonic at 24 Hz is the weakest of these 
lines (not even showing on Fig. 5), and so its proximity to the 25-Hz 
operating frequency has no ill effects. 


4.2.2 Frame currents 


To continue the discussion of noise sources, this section describes 
the problem of distortion of the cable’s field by currents induced in 
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SCARAB’s structure. This type of noise is particularly insidious because 
of its perfect correlation with the signal. Although eddy currents in 
SCARAB’s massive metal objects (such as the manipulator arms and 
flotation tanks) were found to have negligible effect, currents circulat- 
ing in the supporting framework created a severe problem which had 
to be faced. 

SCARAB is constructed as a welded framework of solid aluminum 
rods which support its flotation tanks, thrusters, sensors, and other 
equipment (see Fig. 1). The framework itself can be considered a 
passive electrical circuit which, because of the large cross section and 
high conductivity of the aluminum rods, has an extremely low imped- 
ance. This low impedance has an unfortunate consequence. Time- 
varying magnetic fields, inducing electromotive forces in the frame- 
work in accordance with Faraday’s law, drive sizable currents in the 
frame members. These frame currents in turn cause secondary mag- 
netic fields which, superposed on the primary fields, have the effect of 
distorting the field geometry on which the cable-locating system relies. 

To appreciate the severity of this effect, consider a simplified frame- 
work in the form of a 1-m square composed of solid aluminum rods of 
4-cm diameter (the size of SCARAB’s main frame members). The 
calculated resistance of such a frame is 100 pQ and its inductive 
reactance is 17f uQ (where fis the frequency in Hz). These calculations 
neglect the resistance of the welds and assume a frequency sufficiently 
low so that the skin effect is negligible (not quite true at 25 Hz, but 
nearly so). At frequencies above 6 Hz the framework is primarily 
inductive, at lower frequencies it is primarily resistive. 

Now suppose that an infinite straight wire carries a sinusoidal 
current, situated in the same plane as the square frame but 1 m away 
and parallel to one edge. A straightforward application of Faraday’s 
and Ohm’s laws shows that for frequencies at which the framework’s 
impedance is inductive (i.e., above 6 Hz), the induced frame current is 
about 5 percent of the primary current. Since a magnetometer mounted 
on the framework could easily be 20 times closer to the frame than to 
the wire, the secondary magnetic field seen by the magnetometer could 
easily rival the primary field in magnitude. 

Although SCARAB’s framework is topologically more complicated, 
the same effect is present. To estimate its impact, we ran tests in 
which a current-carrying wire was placed at various locations within 
a meter of the vehicle, and measurements were. made of the magnetic 
fields (primary plus secondary) and of the currents flowing in the 
frame members. The amplitudes of the frame currents, measured with 
a specially modified clamp-on ammeter, were of the same order of 
magnitude as predicted in the preceding paragraph. Consequently, the 
field anomaly was large enough to defeat the cable-locating algorithm. 
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We contemplated three possible solutions to the problem. One 
thought was to reduce the operating frequency, thereby reducing the 
induced EMF proportionately. This solution would require lowering the 
frequency to 1 Hz or below, however, because at frequencies above a 
few hertz the impedance is inductive and therefore decreases along 
with the EMF as the frequency is decreased, so that the currents remain 
the same. Since frequencies at or below 1 Hz are objectionable for 
many reasons (e.g., motional noise), we discarded this alternative. 
Another thought was to mount the magnetometers far away from any 
frame members, thereby reducing the adverse effects of the frame 
currents. In practice, this alternative requires four booms that hold 
the magnetometers one or two meters outside the vehicle’s perimeter 
and retract to prevent damage during vehicle launch or recovery. We 
did not attempt to implement such a system, although it is probably 
feasible, because a simpler alternative was available. 

The final solution was to prevent the frame currents from flowing in 
critical frame members by cutting these members and splicing them 
back together with insulating splices. Obvious choices for splicing were 
the four frame members on which the four magnetometers are 
mounted. Also spliced were the two vertical frame members at the 
front of the vehicle (see Fig. 1). Frame currents flow in the remaining 
frame members, of course, but their distance from the magnetometers 
reduces the adverse effects to acceptable levels. 

The combination of frame currents and thruster noise constrains the 
placement of magnetometers to the front of the vehicle, where there 
are no thrusters and where the density of frame members has been 
minimized to improve TV visibility. Initial plans included the place- 
ment of one magnetometer further aft, but the high density of frame 
members at the intended location (each of which would have had to 
be cut and spliced) forced the relocation of this magnetometer to the 
front of the vehicle. The consequence, as mentioned in a preceding 
section, is that the cable-localization system is blind to a cable which 
lies under the vehicle and perpendicular to its center line, but this 
blind spot does not significantly reduce the usefulness of the system. 


4.2.3 Motional noise 


The strong source of noise evident below 5 Hz in Fig. 5 is motional 
noise, attributable to changes in magnetometer orientation relative to 
the earth’s magnetic field as SCARAB pitches, rolls, and yaws. The 
identification of this noise is verified by its absence when SCARAB is on 
a stable platform (such as the sea bottom). When SCARAB is surfaced 
in a seaway, a spectral peak in the motional noise at the frequency of 
the surface water waves is sometimes identifiable. Above 5 Hz, how- 
ever, motional noise is negligible relative to other noise sources. 
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4.2.4 Broadband noise 


In addition to narrowband noise from the sources described above, 
SCARAB’s noise spectrum includes broadband noise which falls off with 
increasing frequency roughly as 1/f (see Fig. 5). This frequency depen- 
dence matches the Fourier transform of a step function, and indeed if 
one looks at the noise in the time domain (restricting attention to a 
quiet frequency band such as the neighborhood of 25 Hz), one sees a 
series of impulsive transients. This behavior suggests that the broad- 
band noise is due to the combined effects of a number of weak, random, 
impulsive sources. 

One example of such an impulsive source is the actuation of a 
solenoid, of which many are used in SCARAB’s hydraulic system. In 
fact, under quiet circumstances the actuation of SCARAB’s solenoids is 
readily observable by a magnetometer more than 2 m away. Other 
potential impulsive noise sources include any movable equipment with 
residual magnetization, such as SCARAB’s manipulator arms and 35- 
mm camera (whose casing is highly magnetized, causing large broad- 
band transients during panning and tilting). Slight random shifts in 
the position of such equipment would cause broadband noise of the 
kind observed. 

Although the identification of the principal broadband sources re- 
mains speculative, there are a number of sources which can definitely 
be excluded from consideration, since their levels are demonstrably 
insufficient to account for the observed broadband levels. These 
sources include detector noise (generated within the magnetometer 
core), atmospheric noise (from nearby man-made sources or from 
worldwide thunderstorm activity, ducted by the ionosphere), and tape 
noise (on the magnetic tape used for recording test data). 


4.3 Signal-to-noise ratio and cable-localization performance 


With the signal and noise levels known, the performance of the 
cable-locating system can be predicted. Consider a typical scenario in 
which SCARAB is looking for a cable break 100 miles offshore. Figure 4 
shows that the net current at the break is about 40 mA (rms) and Fig. 
5 implies that the noise level in a 1-Hz band at 25 Hz is roughly “y 
(rms). The range at which the cable’s signal, given by eq. (1), equals 
the noise level is 32 m, which may be termed the “cable-detection 
range.” 

At this range, however, and even well inside it, the cable-localization 
algorithm still does not have adequate signal-to-noise ratio to localize 
the cable because it is the field differences (rather than the fields 
themselves) that contain the significant information for localization. 
Therefore, to calculate a “cable-localization range” by equating signal 
and noise, one should use as the signal the field difference AH = 
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(Ap)(0H/dp), where Ap is a measure of the mutual separation of 
SCARAB’s sensors. To allow for the noise contributions (assumed to be 
correlated) from two magnetometers, the noise level should be in- 
creased by a factor of two. With Ap = 1 m, this calculation gives a 
localization range of 4 m. 

AS SCARAB approaches the cable, therefore, it passes through three 
regions in which it can infer increasing levels of information concerning 
the cable. Outside the detection range, it does not sense the cable. 
Inside the detection range but outside the localization range, it detects 
the proximity of the cable but has insufficient information to determine 
the cable’s location. Inside the localization range it can localize the 
cable, with the precision of this localization increasing as the distance 
to the cable decreases. For a net cable current J = 40 mA (rms), the 
detection and localization ranges are 32 m and 4 m, and for other 
values of I they scale as J and I”, respectively. 

Localization of the cable is reasonably precise only at quite modest 
ranges (a few meters or less). The fundamental obstacle to localization 
at greater ranges is SCARAB’S compactness, which imposes severe 
constraints on the size of the magnetometer array. However, the 
description in Section VI of SCARAB’s operational use shows that the 
present localization range is quite adequate. 


V. CABLE-LOCATING SYSTEM: DETAILS OF DATA COLLECTION, 
PROCESSING, AND DISPLAY 


Preceding sections have discussed the cable-locating system in 
rather general terms, omitting details of the actual hardware and 
software which carry out the functions of data collection, processing, 
and display. Some of these specifics will now be discussed, to promote 
a better understanding of how the system works in practice. 


5.1 Data collection 


The sensors are four commercial three-axis ring-core fluxgate mag- 
netometers. Each magnetometer is mounted in a cylindrical aluminum 
pressure housing with three preamplifiers whose function is to remove 
from the magnetometer outputs the dc components (which are pro- 
portional to the three components of the earth’s dc field) and to 
amplify the residuals for transmission to an electronics housing at the 
stern of the vehicle. The relative gain among the three preamplifiers 
is adjusted to compensate for the differential attenuation of the mag- 
netic field components by the cylindrical aluminum container (whose 
eddy currents reduce the 25-Hz field inside to about 82 percent of its 
outside value, with slight attenuation differences between the axial 
component and the other two). 

Each magnetometer signal is then passed through a bandpass filter 
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at 25 Hz and band-elimination filters at 30 Hz and 60 Hz (to eliminate 
the pump-motor and power-frequency noise). The signals are then fed 
to rms detectors whose outputs are sampled five times per second, 
digitized by a 12-bit analog-to-digital converter, and multiplexed onto 
SCARAB’s umbilical cable for transmission to a minicomputer on the 
control ship. 

In addition to the rms amplitudes, the cable-localization algorithm 
requires the relative phases of the three field components from each 
magnetometer to allow determination of the orientation of the vector 
magnetic field. These relative phases (which are either 0° or 180°) are 
continuously determined by phase comparators whose outputs, sam- 
pled five times per second, are digitally coded and packed with the 
digitized amplitudes for transmission to the surface minicomputer. 

The process of digitizing the amplitudes introduces a certain amount 
of error because of the uncertainty in the digitizer’s least-significant 
bit and because of bias and fluctuations in its reference voltage. This 
error normally affects only the two lowest-order bits, and is equivalent 
to magnetic noise of about 0.1 to 0.2 y (rms), somewhat less than (but 
not negligible compared with) the vehicle noise at 25 Hz. This digital 
noise level could be reduced by increasing the preamplifier gain, with 
the penalty that the digitizer, which presently saturates at fields of 
about 120 y (rms), would saturate at proportionately lower levels. The 
present gain setting was predicated on keeping the effects of digital 
noise below those of vehicle noise while maximizing the field at 
saturation. It follows that any future improvement or degradation in 
vehicle noise may require corresponding gain changes to maintain the 
desired relationship between digital noise and vehicle noise. 


5.2 Data processing 


The data collection system just described provides the shipboard 
minicomputer with three field-component amplitudes and their rela- 
tive phases from each of four magnetometers. Sensors on the vehicle 
also provide its digitized heading, pitch, roll, and height above the 
bottom. New data are provided five times per second. 

The first step in the real-time processing of these data, after the 
relative phases are used to assign signs to the rms component ampli- 
tudes, is to calculate a time average of the magnetic fields using an 
exponentially decaying weight function, 


(H) -| Ce H(t), (2) 


where 7 is an averaging time and H(—¢) denotes the magnetic field at 
a time ¢ in the past. This running average, whose purpose is to reduce 
the effects of the noise, is updated with each new input of data (five 
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times per second). The value of 7 is chosen by SCARAB’s operator, and 
may be changed at will as dictated by the operational situation. The 
rationale for particular choices of + will be discussed below. 

Note that this averaging is not the only averaging undergone by the 
data. The 25-Hz bandpass filter (whose Q@ = 20) has the side effect of 
averaging the analog signals according to eq. (2) with rt = 0.25 s. 

From these averaged data the computer determines cable location 
twice per second according to the algorithm described geometrically in 
Section 3.4. To express the algorithm mathematically, let the cable 
location be specified by the vectors ro and a, where Yo is the vector 
from the coordinate origin to the closest point on the cable and ais a 
unit vector parallel to the cable. If the cable carries a net current J, 
then the magnetic field observed at any point r is 


I a X (r — ro) 


H(r) = — ———__—_—__ . 
) 27 |r — Yo — Aa-(r — Yo) |” 


(3) 
This equation is the vector generalization of the scalar eq. (1). 

The core of the data-processing problem is to calculate the cable’s 
location (expressed by ro and a) from two measurements of the 
magnetic field (H,, and H,) made by two magnetometers located at 
coordinates r,, and r,. The dependence of ro and a on the measured 
parameters H,,, H,, rm, and r, is given by the two equations below, 
whose derivation follows: 


= Ym? Hmn H2 ae r,-H, H,,:H, 


= H,, 
ms H?, HH? _ (Hn-H,)? 
r,-H, H?, — Yn: H,,Hn- HH, 
Thy 4 
HH? — (H.-H)? o 
H,, X Ha 

4 = 5 
* TH, X Hal (6) 


These equations can be derived using mathematical arguments that 
parallel the geometrical discussion of Section 3.4. For variety, however, 
we present a more formal derivation. 

Suppose for the moment that H,, and H, are not parallel, in which 
case it is clear that the three vectors H,,, H,, and H,, <x H, comprise 
a basis spanning three-dimensional vector space. The vectors ro and a 
must therefore be expressible as 


Yo = oH» + BHn + YHn X Hn, (6) 
4=AH,, + pH, + *Hn X Ha, (7) 


for some choice of the scalars a, B, y, A, p, v. Since eq. (3) shows that 
a must be perpendicular to both H,, and H,, it follows that A = p = 0. 
Choosing v so that @ is a unit vector results in eq. (5) for a. 
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Turning now to eq. (6) for ro, note that by definition ro is perpendic- 
ular to a (and therefore to H,, X H,,), so y = 0. To find a and £, take 
the inner product of eq. (6) with H,, and H, in turn. The resulting 
scalar equations may be solved for a and B, 


_ Yo*Hn H2 — ro-HaHm-Hn 

0” HA? — (H,-H)? 

Yo-H,, H2, — ro-H,,H.,-H, 
H?,H? — (H,-H,)° 


But eq. (3) shows that (rm — Yo)-Hm = 0, so that ro-Hn = rm-Hn (and 
similarly for n). When this identity is substituted into eqs. (8) and (9), 
which in turn (along with y = 0) are substituted into eq. (6), one 
obtains eq. (4) for ro. 

SCARAB’s computer uses eqs. (4) and (5) as follows. For each pair of 
magnetometers, the computer calculates a value of ro and a, and 
furthermore calculates two estimates of the cable current by twice 
solving eq. (3) for the current necessary to account for each of the two 
observed field magnitudes. Because four magnetometers can be paired 
in six ways, these computations provide six estimates of ro and 4 which 
are then averaged with the weight function, 


|H.,. x H,|4 
Wmn = Ht.Hi 5) 
which expresses the fourth power of the sine of the angle between H,, 


and H, (a choice of Wmn whose motivation was discussed in Section 
3.4). For example, 


(8) 


B= (9) 


(10) 


4 
; » WwyYo(Hi, H,, Yi, r;) 
(10) = i>j=1 


aaa Raa ea (11) 


and similarly for a (which, in addition, is normalized to unit magni- 
tude). An estimated cable current is also obtained, using a similar 
weighted average. In all these averages, magnetometer pairs for which 
H,, and H,, are nearly parallel (to within one degree) are omitted to 
avoid the computational singularity when H,, | | H, [which would cause 
the denominators of eqs. (4) and (5) to vanish]. Also omitted are any 
magnetometers whose contribution SCARAB’s operator wishes to ex- 
clude (e.g., a malfunctioning magnetometer). 

Having thus calculated an estimated cable location, the final com- 
putation uses the heading, pitch, and roll data to transform this cable 
location from a coordinate system fixed in the vehicle to a coordinate 
system independent of vehicle orientation. The cable location is then 
ready to be displayed to SCARAB’s operator. 
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5.3 Data displays 


The computer calculates a new cable-location estimate about twice 
per second, in the manner described above, and communicates this 
estimate to SCARAB’s operator by means of a plasma display panel. 
Both computer and display are located in SCARAB’s operator control 
hut, pictured in Fig. 6. The display format shown in Fig. 7 is one of 
several which the operator can invoke to aid in the various phases of 
SCARAB’s operations. The particular format of Fig. 7 shows the esti- 
mated cable location in plan view relative to SCARAB’s position, with 
SCARAB’s heading indicated both numerically and by the orientation 
of the vehicle symbol (whose rectangular outline is drawn to scale to 
represent the actual location of the vehicle’s perimeter). The estimated 
cable current and the operator-selected averaging time are also dis- 
played. 

When the vehicle is far from the cable, the magnetic data are 
dominated by noise and the cable-location estimates are random and 
meaningless. To suppress such estimates, the computer examines the 
cable-location estimates to establish their credibility before displaying 
them. For example, a cable-location estimate is not displayed if it is 
blatantly inconsistent with the preceding estimate. The criteria for 
credibility are sufficient to identify most estimates that are dominated 
by noise, so that in the absence of a cable the operator is not distracted 
by a random succession of meaningless lines. 

A variation of this display format, shown in Fig. 8, includes a history 
of the estimated burial depth of the cable. This estimate is determined 
by subtracting SCARAB’s height above the bottom (as measured by its 
altimeter) from its estimated height above the cable (given by the 
vertical component of ro, appropriately adjusted if a nonhorizontal a 
indicates a sloping sea bottom). The estimated burial depth is more 
severely affected by noise than the estimated cable location, because 
the former is the difference between two larger quantities, each subject 
to noise. Accordingly, provision has been made for the burial depth to 
be time-averaged, with its averaging time chosen by the operator. This 
averaging time can be chosen substantially longer than the averaging 
time for the raw data, because at slow vehicle speeds the actual burial 
depth changes quite slowly with time. 

A third display format (Fig. 9) shows a history of the three magnetic- 
field components measured by a single magnetometer (selected by the 
operator). This display is useful for troubleshooting, since the effects 
of a malfunctioning magnetometer are immediately apparent. It is also 
useful for cable detection when SCARAB is more than a few meters 
away from the cable (i.e., outside the cable-localization range but 
inside the cable-detection range). Furthermore, in cases when the 
signal is so weak that the cable-localization algorithm is ineffective at 
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Fig. 6—SCARAB operator control hut. Computer keyboard and display are partially obscured but visible at the 


right-hand seat. 
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Fig. 7—Real-time computer display showing estimated cable location relative to 
SCARAB, in plan view. 


any range, cable tracking using this display may still be possible, as 
described in the next section. 


VI. OPERATIONAL USE OF THE SYSTEM 


At the time of this writing, sea trials of the cable-locating system 
have been completed. The first tests were done by emplacing a length 
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Fig. 8—Real-time computer display showing estimated cable location and a short 
history of estimated burial depth. 


THE SCARAB SUBMERSIBLE 561 


SIGNALS FROM MAGNETOMETER NO.1_—s scarap’s 


FWD 4 TRUE 
HEADING 

180 
0 

1-SEC 

DATA 

AFT —4 AVG 
PoRT 4 
0 
STBD -4 
up 4 
0 
DOWN —4 





60 SECONDS ————_———- 


Fig. 9—Real-time computer display showing the three components measured by a 
magnetometer as a function of time. 


of trial cable, applying a 25-Hz tone, and then acquiring and tracking 
the cable using the system. Later, tests were done on two operational 
cable systems. Finally, the system was used during the burial of a 
section of an operational cable. 

Background noise levels were noted under various conditions 
throughout the trials. During vehicle pre-dive checkouts, noise levels 
were noted with the vehicle on deck. In this situation, magnetometer 
outputs were noted to vary randomly in the range of +0.3y. After the 
vehicle was launched it maneuvered on the surface for some time 
without using vertical thrusters and the outputs would decrease to 
+0.2y. At the beginning of a dive, the outputs increased to as much as 
several gammas. These levels result from the 25-Hz field components 
emitted from the vertical thrusters as they pass through the speed at 
which their slip frequency is 25 Hz. These noise levels were evident 
many times during descent as the vertical thruster levels were varied 
by the operators. 

When the vehicle reached the bottom, vertical thrust was increased 
to 100 percent and magnetometer levels subsided to less than +0.5y. 
Transients of about 2y were observed on all channels whenever vehicle 
lights were turned on (or off). 

The location of the dive was generally far enough from the trial 
cable so that acquisition would be initiated by dead reckoning. As the 
vehicle approached the estimated cable position, a numerical display 
of magnetometer signals was observed to monitor the vertical fields 
measured by the magnetometers. (As an alternative, the vertical field 
of a single magnetometer could be observed in a time series using the 
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display of Fig. 8.) When an appreciable vertical field was observed, the 
display of Fig. 6 was used to allow the vehicle operator to track the 
cable. 

When a specially laid trial cable was used an initial tone current of 
200 mA (rms) was used. Field levels as large as 22y were indicated and 
the cable was successfully tracked at an altitude of 4 meters. 

After some confidence was established, tone current was reduced to 
100 mA and then to 50 mA and indications were not noticeably 
degraded. Further reduction of the current to 10 mA (rms) produced 
an objectionable disturbance in the display, but it was still useful. 
There is some evidence that the actual cable currents may have been 
somewhat smaller than these values, because of injection losses where 
the current was applied to the cable; this effect was not quantified. 

The burial-depth feature gave somewhat inaccurate results when 
the vehicle was maneuvering off the bottom, because of fluctuating 
altimeter readouts; however, stable and repeatable results were ob- 
tained when the vehicle was held on the bottom by its thrusters. In 
this mode of operation, the operator notifies the computer that SCARAB 
is being held on the bottom, and the computer therefore ignores the 
altimeter in determining burial depth. Thus, any inaccuracy in the 
altimeter is avoided. In addition, maneuver noise is reduced and 
vertical thruster levels tend to remain above the range in which 25 Hz 
emission occurs. 

During the trials on operational cables, SCARAB was operated about 
100 nmi from the cable terminal where the tone was applied to the 
cable. Tone currents of 20 and 30 mA (rms) were used. As a result, the 
actual 25-Hz current at the point of operation was only about 5 mA 
(rms). This proved marginal. Cable acquisition with the vehicle at an 
altitude of one meter could not be accomplished using the magmeto- 
meter system; however, after the cable was acquired visually, tracking 
was possible. The system indicated an estimated tone of 3-6 mA during 
- these tests. It should be remembered that during acquisition and 
tracking of faulted cables (SCARAB’s principal mission), current levels 
can be used which are an order of magnitude higher than these. 

The system was later used during the burial of an operational cable. 
Here, the indicated tone current was 10 mA and both acquisition and 
tracking were possible. 

In view of the brevity of these trials, a full evaluation of the 
capabilities and limitations of the system must await actual operations. 
However, since the sea trials tended to confirm performance predic- 
tions, it seems likely that the system’s performance is well enough 
known to support the following speculations as to the likely operational 
use of the system during a cable-repair operation. 

With SCARAB on board, the cable-repair ship will steam to the - 
estimated fault location and, perhaps with the aid of towed electrodes 
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or induction coils to sense the 25-Hz tone, will position itself shoreward 
of the fault and relatively near the known cable track (though not so 
near that the direction to the cable is ambiguous). After launch, 
SCARAB will proceed on or close to the seabed at a speed of about 1 
knot in a direction that will cross the cable at right angles, with 
SCARAB’S crew monitoring the computer display of Fig. 9 to watch for 
the cable’s presence. 

At a few tens of meters from the cable, its signal will start to rise out 
of the noise. As the cable is crossed, the characteristic signature shown 
in Fig. 9 will be observed, indicating that the cable has been crossed. 
SCARAB’S operator will then turn the vehicle back toward the cable, 
intercepting it obliquely and then swimming directly above it toward 
the fault, using the display of Fig. 7 to facilitate cable tracking. Except 
in severe Ocean currents, SCARAB should be able to stay within a few 
meters of the cable, so the restricted range of the localization algorithm 
should pose no handicap. Gentle bends in the cable should not affect 
the tracking algorithm adversely, as long as their radius of curvature 
is large compared with the dimensions of the magnetometer array. 
Sharper bends normally will be unburied, therefore causing no confu- 
sion if visibility is adequate. If the cable current is less than about 10 
mA (rms), it may be insufficient for tracking using the tracking 
algorithm. However, SCARAB may still be able to track the cable using 
the display of Fig. 9, either by steering so as to maintain a zero vertical 
field component or by zigzagging across the cable, ascertaining from 
the display each time the cable is crossed. 

When the cable fault is reached, it may very well be recognizable 
visually, either because the cable was pulled out of the bottom when 
it was broken or because the source of the damage (e.g., a trawler’s 
otter boards) left visible tracks in the bottom. If there is no visual 
indication, the fault location may be recognized by the sudden decrease 
in the magnetic field, causing a deterioration of the cable-location 
estimate as the fault location is passed, accompanied by a sudden 
decrease in the estimated cable current. This effect results from the 
nature of the field near the end of the cable, where the field magnitude 
falls to zero in a horizontal distance comparable to the vertical distance 
of the vehicle above the cable. SCARAB’s cable-location estimate will 
therefore deteriorate completely in a distance of about a meter (cor- 
responding to an elapsed time of about 2 s at a speed of 1 knot), 
showing that the fault has been passed. With this realization, SCARAB’sS 
operator can circle back, reacquire the cable, and position SCARAB on 
the bottom near the fault, perhaps recognizing the precise fault location 
by the fact that when SCARAB is directly over the fault the magnetic 
field (and hence the estimated cable current) will have fallen to half 
the value it had before the area of the fault was reached. Once SCARAB 
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is on the bottom at the fault, it will use its dredger to unbury the cable, 
and the job of the cable-localization system is complete. 


Vil. SUMMARY 


Increasing reliance on burial to protect undersea cable has led to 
the construction of the SCARAB submersible, whose ability to find and 
track such cable is clearly crucial to its success. SCARAB’s primary 
cable-finding system is similar to previous ship-based cable-finding 
systems in that it requires a 25-Hz tone on the cable, but it is 
considerably more sophisticated than these previous systems so that 
the cable-tracking job is made as easy as possible. A ship-based 
minicomputer processes signals from four three-axis magnetometers 
aboard the vehicle, estimating the relative cable location, and display- 
ing it pictorially to SCARAB’s operator. The principal difficulties en- 
countered in implementing the system were caused by two types of 
noise source on the vehicle: the thruster motors, and currents induced 
in SCARAB’s framework by the cable’s 25-Hz field. The thruster noise 
was circumvented by placing the magnetometers as far as possible 
from the thrusters and avoiding the thruster speed which would result 
in noise concentrated at 25 Hz. The frame-current problem was re- 
duced to acceptable levels by inserting insulating splices in the frame 
members nearest the magnetometers, after mounting the magnetom- 
eters at locations where frame members were relatively sparse. Once 
these problems were overcome, the system was found to function 
satisfactorily with cable currents down to about 10 mA (rms). So far, 
operational use has been rather limited, and the outer limits of the 
system’s performance have not yet been defined. Nevertheless, it 
seems clear that the system does provide SCARAB with the necessary 
capability for finding and tracking buried cable. 
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APPENDIX 
Low-Frequency Propagation on Undersea Coaxial-Cable Systems 


This appendix concerns the propagation on undersea coaxial-cable 
systems of signals whose frequency is sufficiently low that the electro- 
magnetic fields are not confined to the cable but penetrate substan- 
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tially into the sea water. The purpose of the analysis is to calculate the 
signal levels in the neighborhood of the cable as a function of distance 
offshore, so that the performance of SCARAB’s cable-locating system 
can be predicted. The analysis includes the effects of the repeaters, 
which have an important influence on signal attenuation. 

In practice, an upper limit on the achievable signal level is given by 
the maximum low-frequency current tolerated by the repeaters. This 
maximum is 200 mA (rms) in the most recently installed system (the 
SG system), but is 100 mA (rms) on earlier systems. When in service, 
all systems require that the low-frequency signal be kept below 20 mA 
(rms). These maxima refer to the center-conductor current inserted at 
the shore termination; signal losses in the cable and repeaters govern 
how much of this signal actually reaches the point of interest (i.e., 
SCARAB’s location), and shielding by the cable’s outer conductor and 
armor wires determines how much of this center-conductor current 
can be seen by SCARAB’s sensors. 

The analysis of this problem divides naturally into two stages: a 
waveguide analysis using Maxwell’s equations to determine the prop- 
erties of the propagating modes, and a transmission-line analysis using 
these mode properties to study signal propagation on a repeatered 
cable system. The principal idealizations of the analysis, necessary to 
preserve the cylindrical symmetry of the waveguide problem, are 
twofold. The electromagnetic properties of the sea bottom are taken 
to be identical with those of sea water (insofar as they affect propa- 
gation on the cable), and the effects of armor wires (if present) are 
assumed the same as the effects of a cylindrical shell of equal con- 
ducting cross section. Both assumptions are plausible, but no detailed 
justification of either is attempted. (The theoretical analysis of coaxial 
cables with helically wound wire shields has recently begun to receive 
attention,’ but only in the limit of infinitely thin wires, which is an 
inappropriate limit for the problem at hand because the diameter of 
the armor wires exceeds the skin depth at the frequencies of interest.) 


A.1 WAVEGUIDE ANALYSIS: THE PROPAGATING MODES 


The electromagnetic fields propagating along a cable are governed 
by Maxwell’s equations which, when expressed in cylindrical coordi- 
nates (p, , z), become 


1a 
—-— pH, = €jwkz, (12) 
pap * 
dE. 
—jk.E, - — = —jopH,, (13) 
dp 
jhkeH, = €'jwE,, (14) 
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where the field dependence on time ¢ and longitudinal distance z is 
assumed to be exp(jwt — jkzz), all fields are taken to be independent 
of the azimuthal angle ¢, and 


e’ =€+0/Jw. (15) 


The constitutive parameters are the conductivity o, the permeability 
p, and the permittivity e. 

Equations (12) through (14) relate the three field components H,, 
E,, and E, which comprise cylindrically symmetric TM modes. A 
similar set of equations (not shown) relates the components E,, H,, 
and H, for cylindrically symmetric TE modes. TM and TE modes without 
cylindrical symmetry also exist. However, the cylindrically symmetric 
TM modes are the only propagating modes at the low frequencies of 
interest here, and the other modes are neglected. 

Eliminating the electric fields from eqs. (12) through (14) yields the 
following equation for H;: 


ala a 
(222 p+ bi) H, = 0, (16) 
where 
k, = (w7e/ — kh?) Im(k,) > 0. (17) 


The general solution of eq. (16), 

H, = AJi(k,p) + BYi(Rpp), (18) 
involves the Bessel functions J; and Y; and the arbitrary constants A 
and B. E, can be derived from eqs. (12) and (18), 
kp 


E, = ——[Ado(k,p) + BYo(Rpp) I. (19) 
€Jw 





Coaxial cable consists of a number (N) of homogeneous cylindrical 
layers in which the electromagnetic fields take the form given by eqs. 
(18), (19), and (14) with A, B, k,, and «’ assuming different values in 
each layer. The constants A and B in each of the WN layers are 
constrained to satisfy the 2(N ~ 1) boundary conditions of continuous 
H, and E, at the N — 1 cylindrical interfaces. Furthermore, the 
requirement that the solution be finite at p = 0 and p = & implies that 
B = 0 in the central region and B = JA in the outermost layer. These 
2N constraints imply 2N linear equations that the 2N amplitudes A 
and B must satisfy. To permit a solution, kz must be chosen so that 
the determinant of the 2N x 2N matrix of coefficients vanishes. 

It is easy to show that the number of such eigenvalues kz correspond- 
ing to propagating modes equals the number of insulating layers. If 
the conductors (including the sea water) are perfect, this equality is 
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clear, for electromagnetic waves can propagate in each insulating layer 
independently, each such layer therefore supporting its own single 
propagating TEM mode. In this limit, the coaxial cable comprises a 
number of parallel and independent waveguides isolated from one 
another by perfectly. conducting walls. If the conducting layers have 
finite conductivity, however, the waveguides are no longer indepen- 
dent, since electromagnetic waves can penetrate the conductors, but 
the number of modes remains the same. In this case the modes overlap 
spatially but differ from one another in their relative distribution of 
energy among the various insulating layers as well as in their longitu- 
dinal wave number k,. 

To illustrate the nature of these modes, we modeled a particular 
cable type (sc cable) and solved the eigenvalue problem. Unarmored 
sc cable comprises a central steel strength member (1.06 cm o.d.), a 
copper inner conductor (1.21 cm o.d.), a dielectric (4.32 cm o.d.), a 
copper outer conductor (4.37 cm o.d.), and a protective jacket (5.28 cm 
o.d.). For these calculations, the conductivities of copper, steel, and 
sea water are taken to be 5.9 X 10’, 4.8 x 10°, and 4 mhos/m, 
respectively, the permittivities of the insulators are 2.285€0 (where €o 
is the permittivity of free space), the permittivity of sea water is 81 €, 
and the permeability of steel is 180 po (with po the permeability of free 
space). 

Figure 10 shows the field distributions corresponding to the two 
propagating modes of unarmored sc cable at 25 Hz. Mode 1 is the 
generalization of the TEM mode which, if the conductors were perfect, 
would be confined between the inner and outer conductors, while mode 
2 is the generalization of the TEM mode which would be confined 
between the outer conductor and the sea. The attenuation lengths 
(defined as —1/Im[k-.]) of these modes are 330 nautical miles and 115 
nautical miles at 25 Hz. At neighboring frequencies, the attenuation 
lengths scale roughly as the inverse square root of the frequency. 


A.2 TRANSMISSION-LINE ANALYSIS 


To predict signal propagation on a cable with N conductors (where 
N = 2 for unarmored cable and N = 3 for single-armored cable), we 
must calculate the appropriate amplitudes of the 2N modes (N types 
of modes, two directions of propagation) on each repeater section. The 
superposition of these modes must satisfy the relevant boundary 
conditions (expressed as constraints on voltages and currents) at each 
end of the section. This analysis requires an equivalent circuit for the 
repeaters, which is shown in Fig. 11.° The repeater’s equivalent imped- 
ance depends strongly on whether dc power is applied (as indicated in 
Fig. 11). The cable’s outer conductor and armor wires (if present) are 
grounded to sea at each repeater, as shown. 
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Fig. 10—Field amplitudes of the two propagation modes of unarmored sc cable at 25 
Hz. The broken lines indicate the boundaries of the cable’s layers: central strength 
member, inner conductor, dielectric, outer conductor, outer jacket, and sea water. 


The transmission-line analysis is largely an exercise in matrix ma- 
nipulation, and its details (which involve considerable notational com- 
plexity) will not be given in the interests of brevity. The procedure 
starts at the offshore cable break, where all N conductors are grounded 
to sea. These N constraints on the mode amplitudes, along with the N 
— 1 constraints that all conductors except the center conductor are 
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Vig. 11—Equivalent circuit of sc repeater at low frequencies. The equivalent resist- 
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THE SCARAB SUBMERSIBLE 569 


grounded to sea at the offshore end of the most offshore repeater, 
provide 2N — 1 constraints on the 2N modes propagating in the last 
cable section (i.e., the one between the last repeater and the cable 
break). The relative amplitudes of these modes can therefore be 
determined, leaving only an arbitrary multiplicative constant (common 
to all the modes) to be adjusted later. These mode amplitudes can be 
manipulated to provide the voltage and current at the offshore port of 
the last repeater, and conventional circuit analysis then provides the 
voltage and current at the shoreward port. These values of the voltage 
and current provide two constraints on the 2N modes in the next 
shoreward cable section, with the other 2N — 2 constraints coming 
from the fact that all conductors but one (the center conductor) are 
grounded to sea at both ends of the cable section. The 2N mode 
amplitudes in the next-to-last cable section can therefore be found. 

In this way, the solution is propagated from the cable break back 
toward the shore terminal, with the mode amplitudes being determined 
in each cable section so that the appropriate constraints on voltages 
and currents are satisfied at the cable break and at both ports of every 
repeater. Once the solution is complete, it is adjusted by a multiplica- 
tive constant so as to satisfy the driving-point condition (e.g., that the 
center-conductor current be 200 mA at the shore terminal). 

Figure 12 shows the results of such a calculation, and assumes a 
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Fig. 12—Propagation of low-frequency signals on unarmored sc cable with a cable 


break 100 nmi offshore. The center-conductor current is 200 mA at the shore point, for 
all frequencies. The step discontinuities occur at the repeaters, which are unpowered. 
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Fig. 13—Propagation of 25-Hz signals on unarmored SG cable, with cable breaks at 
various locations. The unlabeled curve extending to 220 nmi is for a cable with no break. 
The center-conductor current is 200 mA at the shore point. The step discontinuities 
occur at the repeaters, which are unpowered. 


cable break to be 100 miles offshore. The repeaters are spaced at 5- 
mile intervals, and their presence is evident by the step discontinuities 
in the net cable current, which is defined as the algebraic sum of all 
currents on the cable’s conductors. The step discontinuities are caused 
by the exchange of current between the outer conductor and the sea 
at each repeater. The center-conductor current itself (not shown) 
decays smoothly with distance, without appreciable step discontinui- 
ties, because the repeaters provide no low-impedance path from the 
center conductor to sea. The center-conductor current is 200 mA at 
the shore point in Fig. 12, but the net current is somewhat less because 
of the return current carried by the outer conductor. This outer- 
conductor current increases with frequency, so that the net current 
decreases with frequency. 

Just as Fig. 12 illustrates the frequency dependence of the signal 
propagation, Fig. 13 illustrates the dependence on the distance to the 
cable break at SCARAB’s chosen frequency of 25 Hz. One predominant 
feature of Fig. 13 is that the cable break reflects the propagating 
modes, doubling the signal amplitude relative to its value in the 
absence of a break. Another important conclusion supported by Fig. 
13 is that the signal losses are primarily introduced by the repeaters, 
rather than by losses in the cable. This conclusion is borne out by the 
fact that the decay from 165 mA at the shore point to 1.7 mA at 220 
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nmi corresponds to an attenuation length of about 50 nmi, much less 
than the 330-nmi attenuation length of the cable’s primary mode. 
Since the repeaters account for most of the signal attenuation (at least 
at 25 Hz), it follows that an improvement in propagation would result 
from a reduction in the repeater’s impedance at 25 Hz. Because 
powering the repeaters reduces their impedance considerably (as Fig. 
11 points out), 25-Hz propagation on a powered system (illustrated in 
Fig. 14) is much less attenuated than on an unpowered system. For 
example, the signal at 200 nmi is increased by more than a factor of 10 
when the system is powered, as comparison of Figs. 13 and 14 reveals. 
Substantial further improvement in 25-Hz propagation on future cable 
systems would result from the inclusion of a bypass network allowing 
the 25-Hz signal to pass through the repeaters unattenuated. Such an 
improvement would increase the usefulness of SCARAB’s cable-locali- 
zation system by allowing useful signal levels to propagate farther 
offshore. 

The calculations in this appendix so far have dealt with unarmored 
cable. Since many existing cable systems (and probably all future cable 
systems) employ armored cable on the continental shelves, it is of 
interest to see how armored and unarmored cable differ in their 
propagation characteristics. To this end, calculations involving ar- 
mored cable have been carried out, modeling single-armored sc cable 
as identical to unarmored sc cable with the addition of a steel armor 
layer (6.3 cm o.d.) sandwiched between two jacket layers (5.75 and 6.76 
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Fig. 14—Same as Fig. 13 except that the repeaters are powered, which reduces their 
impedance at 25 Hz. 
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cm o.d.). The three propagation modes of armored cable have atten- 
uation distances at 25 Hz of 260, 140, and 86 nmi (compared with 330 
and 115 nmi for unarmored cable). Propagation on armored cable at 
various frequencies is illustrated in Fig. 3, and propagation at 25 Hz 
with cable breaks at various locations is illustrated in Fig. 4. Compar- 
ison of Figs. 3 and 4 with their counterparts for unarmored cable (Figs. 
12 and 13) shows one outstanding difference: the net current on the 
armored cable is reduced by the shielding effect of the armor layer 
(which carries a substantial amount of return current, particularly at 
the higher frequencies). At 25 Hz, for example, the net current on 
armored cable is only about 40 percent of that on unarmored cable. 
This ratio is about equal to the quantity exp(—Ap/65), where Ap is the 
thickness of the armor layer and 6 = (2/wyc)'” is the skin depth in 
steel (which is a better shield even than copper because of its high 
permeability », which more than compensates for its lower conductiv- 
ity o). The attenuation of the signal with distance offshore, being 
governed primarily by the repeaters rather than the cable, is the same 
for both armored and unarmored cable. 


A.3 SUMMARY 


This appendix has presented the results of theoretical calculations 
of propagation on sG repeatered cable systems (both armored and 
unarmored). Propagation at various frequencies, with cable breaks at 
various locations, and with both powered and unpowered repeaters, is 
displayed in Figs. 3, 4, 12, 18, and 14. These figures show that, if 10 mA 
(rms) is taken as the necessary current level for cable localization by 
SCARAB, and if 200 mA (rms) of 25-Hz center-conductor current is 
applied at the shore terminal of an unpowered unarmored SG system, 
then breaks as far offshore as 170 nmi can be found (with this value 
reduced to 125 nmi for armored cable, because of the shielding effect 
of the armor wires). These ranges can be increased if dc power can be 
applied to the repeaters, thereby reducing their impedance. If less 
current is applied at the shore point (as on earlier cable systems, which 
can tolerate only 100 mA, or on operational systems, which can tolerate 
only 20 mA), the maximum effective range of SCARAB’s cable-locating 
system is correspondingly reduced. 

These operational ranges could be increased for future cable systems 
by designing the repeaters with a separate low-impedance path for the 
25-Hz signal to avoid the constraints and losses imposed on this signal 
as it traverses the dc power path. 
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